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Abstract. We analyze the long time behavior of solutions of the Schrodingcr 
equation iipt = (—A — b/r + V(t,x))ip, x £ R 3 , r = \x\, describing a Coulomb 
system subjected to a spatially compactly supported time periodic potential 
V(t, x) = Vit + 2ir/u>, x) with zero time average. 

We show that, for any V(t,x) of the form 2Q(r)s'm(ujt — 9), with H(r) 
nonzero on its support, Floquet bound states do not exist. This implies that 
the system ionizes, i.e. P(t,K) = J K \ip(t, x)\ 2 dx — ¥ as t — > oo for any 
compact set i(Cl 3 . Furthermore, if the initial state is compactly supported 
and has only finitely many spherical harmonic modes, then P(t, K) decays like 
t~ 5 / 3 as t — > oo. 

To prove these statements, we develop a rigorous WKB theory for infinite 
systems of ordinary differential equations. 

Keywords: Ionization, time periodic Schrodinger equation, resonances, Flo- 
quet theory, Coulomb potential. 
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1. Introduction and overview of results 

The long time behavior of solutions of the Schrodinger equation of a system with 
both discrete and continuous spectrum subjected to a time periodic potential is a 
longstanding problem. Powerful results have been obtained under various assump- 
tions on the potentials, see [HI El LTl [HI HI1 1321 GUI (Ml EI]> and references therein. 
In particular, there are conditional results on the ionization of the Hydrogen atom, 
subjected to an external time-harmonic dipole field V(t, x) — E ■ x cos tut if E is 
sufficiently small, see (43] [44] . In addition, Moller and Skibsted proved the equiva- 
lence of absence of point spectrum and ionization for a large class of such systems 
subject to periodic fields [32]. There are also detailed results about the behavior of 
the wave function for systems subjected to general time periodic potentials, decay- 
ing faster than r~ 2 , under the additional assumption of absence of point spectrum 
of the Floquet operator, see [20] . 

None of these results however prove or disprove ionization of Coulomb- bound 
particles subject to time-periodic forcing of fixed amplitude and zero average. In 
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fact, such results have only recently been obtained even for simple model systems, 
see [TTJ [T^l [T31 [THl [HI I3U] and references cited there. For a periodic dipole field 
of nonzero average ionization was proved in |33| (we note that the time averaged 
Hamiltonian has no bound states in this seetting). 

What experiments and simplified models show is that the behavior of systems 
with both discrete and continuous spectrum, subject to time-periodic fields of ar- 
bitrary strength, can be very complicated. For amplitudes where perturbation 
theory is not applicable (such fields are becoming of increasing practical impor- 
tance in technology), qualitative departures from the behavior at small fields are 
observed. There are even situations, see e.g. |12) . where for small enough fields 
ionization occurs for all initial states while for larger fields there exist localized 
time- quasiperiodic solutions of the Schrodinger equation, i.e. Floquet bound states. 
Though these situations are rather exceptional, constructive methods of analysis 
are required to determine the outcome in specific settings. 

In this paper we prove ionization for Coulomb systems with very special (non- 
dipole) type of forcings of arbitrary magnitude. This is equivalent to establishing 
the absence of point and singular continuous spectrum of the corresponding Floquet 
operators. We also obtain the large time behavior of the wave function. The 
time decay of the wave function, for compactly supported initial conditions, is of 
order i~ 5 / 6 . This differs from the i~ 3 / 2 or, exceptionally, i -1 / 2 power law found 
for shorter range reference potentials, see [20j[T5]. The nonperturbative methods 
include the development of rigorous WKB techniques for infinite systems of ODEs. 

1.1. The Coulomb Hamiltonian. In units such that h 2 /2m = 1, the Coulomb 
quantum Hamiltonian of a Hydrogen atom (more generally a Rydberg atom) is 

(1) He = -A - - 

r 

where 6>0, r — \x\, x ^ M. 3 and A is the Laplacian. It is well known, see e.g. 
|28j . that He is self-adjoint on the Sobolev space H 2 (M. 3 ) = D(— A), the domain 
of —A (cf. also [5S], p. 303). The spectrum of He consists of isolated eigenvalues 
E n = —b 2 /4n 2 , with multiplicity n 2 , and an absolutely continuous part, [0, oo). 

1.2. Setting. Our starting point is the time evolution of the wave function ip(t, x) 
of the Hydrogen atom described by the Schrodinger equation 

(2) iip t = H c iP + V(t,x)il;; V(0, x) = M*) e H 2 {R 3 ) 

where V(t, x) = ^ ilj(x)e l: ' u}t is real valued and = 

The operator He + V(t,x) satisfies the assumptions of Theorem X.71, p. 290, 
in [3T] v.2.; Theorem X.70, p. 285 also applies in our setting. Thus, for any t, 
ip(t, •) € i? 2 (R 3 ), and the unitary propagator U(t) for @ is strongly differentiable 
in t; see £15.11 for a short proof in our case. 

Assumption 1. The Q,j(x), j € Z are smooth inside a common compact support, 
chosen without loss of generality to be the ball Bi C M 3 of radius 1, and + 

liDII^IU-tB!) < oo. 
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1.3. Ionization. We say that the system ionizes if the probability to find the 
particle in any compact set vanishes for large t, i.e., for any a > we have 

(3) P(t, B a ) = / \ip(t,x)\ 2 dx as t ->• oo 

where B a = {x : |x| < a}. To prove ionization, it clearly suffices to prove Q for all 
a > 1. 

A simple way in which ionization may fail is the existence of a solution of the 
Schrodinger equation in the form 

(4) ip(t, x) = e l<pt v(t, x) with 4> e R and v G L 2 ([0, 2tt/lu] x M 3 ) time-periodic. 
Substitution in ([2]) leads to the equation: 

(5) Kv = 4>v 
where 

(6) K = i^-(-A-br- 1 + V(t,x)) 

is the Floquet operator, densely defined on L 2 ([0,2tt/lu] xM 3 ); 0^»el 2 implies 
by definition that <fi £ <r p (K), the point spectrum of K. 

Somewhat surprisingly, in all studied systems, o~ p {K) ^ is in fact the only 
possibility for ionization to fail. As we will show this is also true for @. The proof 
of ionization also implies that K does not have any singular continuous spectrum. 
This turns out to be a consequence of the existence of an underlying compact 
operator formulation, the operator being closely related to K. Generic ionization 
is then expected since L 2 solutions of the Schrodinger equation of the special form 
dU are unlikely. We prove that for V(t, x) = 2fi(r) sin(W — 6), VL > on [0, 1] and 
sufficiently smooth, they do not exist. 

1.4. Laplace space formulation. For ijj g H 2 (M. 3 ), the Laplace transform 

/•OO 

$(p, •) := / i/j(t, -)e- pt dt 



exists for p£l, the right half complex plane, and the map p — > ijj(p, •) is H 2 valued 
analytic in Re p > 0. The Laplace transform converts the asymptotic problem (J3]) 
into an analytical one. 

To improve the decay in p of the Laplace transform, it is convenient to write 

(7) i/>(t,x) =^ (x)e- t + y{t,x) 

Now, y(t,x) satisfies 

(8) 

iyt - HcV - V(t, x)y = e" 4 + H c ^o + V(t, x)^ ] = -y°(t, x); y(Q, x) = 

Standard arguments (see Appendix 15 .2p show that the t~ Laplace transform of y, y 
is in H 2 and satisfies 

(9) (H c - ip)y(p, x) = y°(p, ^j(x)y(p - ijw, x) 

where 



(10) y°(p, x) = -— + HciPo - J2 J 



l+p\ ' j ^—<l+p-ijoj 
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1.5. The homogeneous equation and the PDE-difference equations. The 

homogeneous system associated to © is 

(11) (— A — b/r — ip)w(p, x) = — Qj(x)w(p — ijuj, x) 

Note 2. (i) Clearly, © and (fTTjl couple two values of p only if (pi — P2) € iwl, 
and are effectively infinite systems of partial differential equations. Setting 

(12) p = p\ + into, with p\ € Cmod (iu>) 
we denote 

y n (pi,x) = y(pi +inuj,x),y°(pi,x) = y°(pi + inu), x), w n (pi, x) = w{pi + inu),x) 
Eqns. Q and (fTTj) now become 

(13) (H c - ipi + nuj)y n = y° - Qj(x)y n -j 

jez 

(14) (He - ipi + noj)w n = - Slj(x)w n -j 

jez 

Note 3. Seen as a differential difference equation, the solution y(p, x) is then a 
vector {y n (pi,x)}„ e z and the whole problem depends only parametrically on p\. 
We have 

(15) y n (p!+iuj,x) =y n+1 (p 1 ,x) 
and the analysis can be restricted to 

§ = {pel: Imp £ [0,cj)} 

where H is the closure of EL There is arbitrariness in the choice of So and, to see 
analyticity in p\ on <9So, it is convenient to allow p\ G H, using (| 1 5|) to identify 
different strips of width ui. 

2. Main results 

Theorem 1. Assume V(t, x) = 2Q,(r) sm(ujt-9), with fl(r) = for r > 1, Q(r) > 
for r < 1 and ^i(r) € C°°[0, 1]. Then a p (K) — and ionization always occurs. 
Furthermore, if ipo(x) is compactly supported and has only finitely many spherical 
harmonics, then P(t, B a ) = 0(t~ 5 / 3 ). 

For the proof, given in §3.10[ £13.111 and M.6\ we develop a relatively general 
rigorous WKB theory for infinite systems of differential equations. This yields the 
asymptotic behavior of w n as n — > — oo. The argument relies on Theorems [5] and [3] 
below. 

Remark 4. The condition 0(1~) ^ simplifies the arguments but these could ac- 
commodate an algebraically vanishing f2. (We also note that some one-dimensional 
models with rough such as a S mass show failure of ionization, see Q2] , [30] and 

m\) 

We will later derive equivalent systems of integral equations, (|2"2"|) . allowing for 
a compact operator reformulation of the problem. 
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Theorem 2. In the setting assuming spherical symmetry in x of the forcing 

V(t, x), ionization occurs iff for allpi € H, has only zero H 2 solutions decaying 
in r v^A . This is true iff o~ p (K) = 0. 

This extends results about absence of singular continuous spectrum of K, [20] . 
to this class of systems, with Coulombic potential and nonanalytic forcing. 
The proof is given in i j3.2l and i j3.8l 

Properties of Floquet bound states for general compactly supported 

V(t,x). 

Theorem 3. If there exists an H 2 nonzero solution w of {1$ decaying in n, 21, 
then it has the further property 

(16) w n — XBxWn for all n < 

with Xa the characteristic function of the set A. 

The general idea of the proof is explained in § 13.91 and the details are given in § 14.71 

Note 5. (i) The Sobolev embedding theorem implies that w n is continuous in x. 
From (fT"4l . w n is piecewise C 2 , implying continuity of Vro„ up to dBi. 
(ii) Equation (1161) makes the second order system (|14[) formally overdetermined 
since the regularity of w in x imposes both Dirichlet and Neumann conditions on 
9Bi for n < 0. Nontrivial solutions are not, in general, expected to exist. 

3. Proofs 

Outline of the ideas. As in our previous work PUJ-Ql)], summarized in [TB] on 
simpler systems, we rely on a modified Frcdholm theory to prove a dichotomy: there 
are bound Floquet states, or the system gets ionized. Mathematically the Coulomb 
potential introduces a number of substantial difficulties compared to the potentials 
considered before (for references, see e.g. [E]), due to its singular behavior at the 
origin and, more importantly, its very slow decay at infinity. 

The slow decay translates into potential-specific corrections at infinity, and stan- 
dard general methods to show compactness in weighted spaces of the Floquet re- 
solvent, such as those in [2U] or our previous ones do not apply. Instead, the 
asymptotic behavior in the far field of the resolvent has to be calculated in detail. 
The accumulation of eigenvalues of increasing multiplicity at the top of the discrete 
spectrum of He produces an essential singularity at zero of the Floquet resolvent 
with a local expansion of the form l p'/ 2 e «M(- J p) 1 f 01 - small p when Rep > 0, 
where A = nb/2. For sufficiently rapidly decaying potentials the exponentials would 
be absent. Their presence clearly makes the analysis at p — of the Floquet resol- 
vent more delicate and is responsible for the change in the large time asymptotic 
behavior of the wave function, from i~ 3 / 2 to i~ 5 / 6 . 

We introduce an extended parameter X — (p 1 / 2 , e lA ^~ lp ^ 1 ) and prove analyt- 
icity of the solution y in X, whose p-counterpart is p small, Rep > 0, and similarly 
in regions near the special points p £ iufL. We reformulate the problem in terms 
of an integral operator £, defined in H3.41 closely related to the Floquet resolvent, 
shown to be compact in a suitable space and analytic in a variable corresponding 
to X. Then, by the Fredholm alternative, (I — £) _1 is meromorphic, and in fact 
analytic in X, since we show absence of eigenvalues of I — £ for any pGi. 



WFor precise conditions, see i|3.4l below and the integral form (l-'2t . 
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3.1. The Hilbert space %. Let % be the Hilbert space of sequences Y = {y„}„ e z, y n & 
L 2 (B a ), with a > 1, and with 

Hill 2 := \\Y\\l = + M) 4/3 |lifell£»(B.) < oo 

Note 6. TTie properties of (I — as Repi — > 0+ ensure that Y(pi, ■) <E HU H 2 
and is locally integrable in p\ along iR. 

We then extend the stationary phase method to such a setting, cf. ^5.71 to 
evaluate, asymptotically for large t, the inverse Laplace transform of y on iM. and 
obtain the ionization result and time decay estimates. 

To show ionization we then have to rule out the existence of a point spectrum of 
the Floquet operator, that is the existence of nontrivial solutions of (HU). We use 
the general criterion in Theorem [3] to show that, if there exists a nonzero solution to 
(|14[) . then a subsequence of {w n }„ e _N would be singular at x = 0, in contradiction 
with NoteEJ (i). 

To find the behavior of solutions for large n, we develop a WKB theory for infinite 
systems of ODEs and find the asymptotic behavior of w n in n in detail. The formal 
WKB calculation of the behavior is straightforward algebra, relatively easy even in 
much more general settings, see A5.9I Justifying the procedure is however delicate, 
and a good part of the paper is devoted to that; cf. §4. lit §4.121 

The procedure of introducing an enlarged set of parameters with respect to which 
the solution is regular, when this does not hold in the original parameter, should 
also be applicable to other problems where complicated singularities arise. 

3.2. Proof of Theorem [2j We show that y has a limit in L] oc on dW = iR where 
it is smooth except for possible poles and a discrete set of essential (but L 1 ) singu- 
larities. Poles arc present iff the integral form (|2"2"|) of (fTl} has nontrivial solutions 
in %. There is sufficient decay in p at infinity, so that, when poles are absent, 
the Riemann-Lebesgue lemma applies, implying that y decays as t — > oo, prov- 
ing ionization -since ipo(x)e~ t obviously goes to zero in this limit. More detailed 
analysis of the resolvent reveals the nature of the essential singularity at p = iufL. 
Stationary phase analysis shows a i~ 5 / 6 decay of the wave function if the initial 
condition is spatially compactly supported and contains only a finite number of 
spherical harmonics. 

Proposition 7. Ionization holds for every i/jq € L 2 iff it holds for any ipo in a set 
densely spanning L 2 . 

Proof. We make use of the standard triangle inequality argument to estimate 
U(t)i/jo, where U(t) is the unitary operator associated to the Schrodinger evolu- 
tion ©. | 

We choose ip in a dense set C^°(K. 3 ), the smooth, compactly supported functions 
in K 3 . Define as usual the angular momentum operators 

2 _ d 2 1 d 1 d 2 
~d9 2+ t^OdO + sin 2 9 W 

and 
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Let Pi :m be the orthogonal projector on {<fi : ~L 2 <f> = 1(1 + 1)0, L z tf> = mcf>} for some 
to e Z, |to| < Z € Nu {0}. 

Since Pi, m = /, we can now assume without loss of generality that tpo € 

P;, m (c c °°(]R 3 )) if /and to are arbitrary. Likewise, if P(t, B a ) decays like t 5 / 3 when 

V->o G -R,m fC£° (M 3 )) , then the same decay rate clearly holds for any tpa given by a 

finite linear combination over (l,m) (but not, in general, for any ipQ € L 2 (R 3 )). 
Further notations. An index of commonly used notations together with their 
pages of definition is given at the end of the paper. As usual we write D e = {z : 
\z\ < e}, B = Di and we denote D+ = D e n {z : argz <E (-tt/4, tt/4)}. We also let 
l t = i[-e,0], H+ c = {p + c : p e H}, i a = {p : Re(p) > 0,lm(p) = a} and for a 
set A, Ay a = A \ £ a . We denote V = H U iIR + , 0(D) will denote some small open 
neighborhood of T>. 

3.3. Step 1. Compact operator reformulation. To investigate the analytic 
properties of ip it is convenient to introduce a new operator Ap which is a complex 
perturbation of He, having no real eigenvalues. More precisely, define 

(17) Ap ■= He — if3(p)XB a (r) — ip; (with the understanding that Ao = Hc — ip) 

where a > 1 and 



(18) /3 = (3(p) 



c > if Imp 6 [— e,p c ] and Rep > 
otherwise 



Here e < uj/2 is small as required in Proposition [171 below, and we choose p c so 
that Pc/uj ^ Z and p c > — _Eo = o 2 /4, the ground state energy of the unperturbed 
atom. Clearly Ap is defined on D(H ) and A* 3 = A-p + ip* + ip. We rewrite ((Tc 
and (|14p in the equivalent form 



(19) [£T C - i Pl +nuj- i/3(p)X Ba {r)} y n = y ° - iP(p)X Ba {r)Vn - J2 n j(x)y n -j 

(20) [H c - ipi +nuj- i(3(p)X Ba (r)} w n = -i(3(p)X Ba (r)w n -y^Clj(x)w n -j 

We show next that A^ 1 is analytic in p € H \ {£ Pc U ^_ e }, and sufficiently regular 
on ?R. Since the parameter p c is artificial, the non-analyticity at £ Pc U of ^l^" 1 
is not reflected in the actual solution y, as discussed in Note [3 

Proposition 8. There exists an open neighborhood O of "D\ {£ Pc U £- e }, not con- 
taining the origin 0, such that the operator 9lp = A^ 1 exists and is analytic 
in p £ O \ (£ Pc U £_ e )- Furthermore, for any p for which 9\p exists, we have 
9%3 : L 2 (R 3 ) -> H 2 (R 3 ). 

The proof is given in i j4.ll 

3.4. Restriction to a ball B; Definition of £. To study ionization, we only need 
to know y(t, x) for X in a fixed (but arbitrary) ball B a D Bi. Henceforth, to simplify 
the notation, we write B a = B. We shall therefore need to study the properties 
of X B ^\pX B . This sandwiched operator (which preserves information about L 2 (R 3 ) 
through built-in boundary conditions on <9B) is the one that we shall most often 
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use below. We recall that p — p\ + inu and y(pi + inu3,x) — y n (pi,x). Since 
flj(x) = XB^ij(x), (TLB")) implies that for igB, 

(21) y n (j>i,x) = XB^isVn + ^b%X b - i/3y n (pi ,x)-J~]toj (x)y n -j (pi , x) , 

jGZ 

where we may assume that B contains the support of i/}q(x), and therefore of 
Note that 9^ depends on n through p = inu + p\. Corresponding to (|2lj). we 
obtain the homogeneous system: 

(22) iy„(pi,x) = Xb^/sXb - i/3w n (pi, x)-}] Clj(x)w n -j(pi,x) 

The elements of % will be denoted by capital letters e.g. {y n } n£Z =■ Y, 
{Vn} n pz ~ : ^ e define the operators T on L 2 (B) by 

{«o}„ = ^B%y° 

and £ on H by 

{£^}„ = Xb^hXb - iPy n (Pi, x) - ^ nj-(g)j/ w _j-(pi, x) 

jGZ 

Then, we rewrite (1211) in the form 

(23) y = XY + £y 

Note 9. We shall see that for any f3 satisfying (|18|). eq. (f2"3"|) has a unique solution in 
H (call it now yW). Thus, away from the artificial cuts, all these solutions coincide 
(since the domain of Y^ corresponding to c = 0, contains all of the others). Hence, 
wherever some has analytic continuation, so will Y^°\ 



The homogeneous system corresponding to (|23p is given by 
(24) w = €w 



Note 10. We have shown (cf. §l-4\ and Q5.ty) that tp and the Laplace transform 
v(Pt') = £ — V*oe ) exist for Rep > 0. The corresponding Y = {y n }, Vn = 
y(inuj +pi, ■), restricted to B, will therefore satisfy i2fy) for j3 — when Repi > 0. 
It will be shown that (S3\) has a unique solution Y € % for any Repi > 0, and that 
Y is analytic in pi £ H and has an Lj oc limit on zR, with sufficient decay in n. 
The implied decay and regularity properties of yip, ■) on iR show that £~ 1 y + e~ t ipo 
(the integration contour taken to be iR) equals ip for x £ B. 

Proposition 11 (Asymptotic behavior of XbJH/jXb)- IfRcp — and \lmp\ —> oo 
(see Note\$), then ||X B %X B || = 0{\p\-^ 2 ) (recall that /3 = if |Im p\ is large). 
Moreover, for any e > 0, Xb^/sXb is analytic in p in an open set containing 
—i(e, oo). 

For Im p — > +oo, the |p| _1//2 decay rate follows from the spectral theorem since we 
are outside the spectrum, while for Im p — > — oo, the rate is obtained using Mourre 
estimates [27] , Theorem 6.1. The rest of the proof involving analyticity is given in 
£14.31 and relies on an explicit representation of the resolvent for He, see £)5.4l Using 
spherical symmetry, the explicit Green's function could be avoided, but in view of 
possible future generalizations to non-spherical V(t, x), we prefer this more delicate 
approach. 
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Lemma 12. We have 1Y € %. The operators & := Y — > {J2j<=z ty( x )Vn-j(pi, x)} ne z 
and £ are bounded in %. 

Proof. We note from PropositionHTIthat — O (|p| -1 ^ 2 ) for large p, i.e. O (In^ 1 / 2 ) 
for large |n|, since p — into +p%. Therefore, from the expression of y° in (I10p . 



^(i + H) 4 / 3 ||{Tr }j| 2 

<^(l + |n|) 4 / 3 (||XB%X B | 



|^o||l» + ||A^o|| 



\x 



B%x B ||||VolU^^ 



|1 +Pi + inuj\ 
l%IU 2 



< C 



\l+Pi +i(n- j)ui\ 



i + Ed + H) 6/7 (ett 

n£Z \jSZ ^ 



|%|| L 2 



il) 



Using (7.12) and (7.13) in [T5] with 7 = 6/7, the above is finite since 

+ bi) 3 / 7 ii^n L2 < 53(1 + \j\m 3 \\ L , < 00 

jez jez 

The proof that © is the same as that of Lemma 27 of [15] , with 7 = | , replacing 
absolute values by norms in x. Since is uniformly bounded (in the operator 
norm) and acts diagonally in n, £ is bounded too. | 



Lemma 13. Both TYo and the operator £ are analytic in pi for 

PiGO(I\ {{£ Pc + iujl} U {£- e + iujZ} U {1 £ + iuZ})) 

Proof. Propositions IBl and ITTI imply that 91^ is analytic in p g O \ {£ Pc U i?_ e } and 
in an open set containing — i(e, 00). Analyticity of 1Y and £ follow from their 
definition (we note £ is a norm limit of analytic operators: its restrictions to the 
subspaccs with nonzero components for |n| < N only). | 



Remark 14. As shown later, (I — £) is invertible. Since the solution Y cannot 
depend on the arbitrary parameters e and p c (see Note\$j), the non- analyticity of £ 
and IYq for 

Pi e {£ Pc + iujZ} U {£- t + iojZ} U {l t + iujZ} 
is not reflected in Y . 

Proposition 15. For Repi > large enough, 123\) has a unique solution in H. 
The inverse Laplace transform in p of y(p,x) =: y n {pi,x), where p — inu) + p\, 
solves the initial value problem |3|) in B (see Note \10\) . 

The proof is given in ij4.2l 
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3.5. Step 2. Regularity of 9Kp t i, m at p = and of £; >m at p\ = 0. 
Define Vtp,i, m = -P;, m % and £;, m = Pjj, m £. 

Note 16 (Compactness versus regularity of Ot^j m ). The term — i/3Xs was intro- 
duced in ^3.21 to ensure that $H^,/, m is bounded in H. Since — i/3Xb is localized in 
x, the shifts in the poles created by the point spectrum of He are smaller as p — > 
(the size of the orbitals of the Hydrogen atom grows when the energy approaches 
zero.) The resulting integral operators have an essential singularity at p = 0. The 
factor Xb is needed to ensure compactness, simplifying the analysis. 

The poles of the resolvent 9l/3 t i >m accumulate at p = from — H, along a curve 
tangent to the positive imaginary p-axis (see Note 1581 in H5.6p . As a result, while 
being uniformly bounded, %Kfj t i t m is not continuous along the imaginary p line at 
zero but oscillates without limit. Boundedness of XbSH^i^Xb (which is not difficult 
to prove) does not ensure boundedness of the solution Y. However, we do have 
analyticity in an extended, two-dimensional, parameter. Let A := \J—vp (with the 
usual branch of the square root, ImA < if p € H) and let X := (p 1 ^ 2 , Z) with 

(25) Z = e^ 

(The dependence of Z on A reflects the actual behavior of the solution.) The 
resolvent is analytic in X and a useful Fredholm alternative can be applied. 

For any a > 1 we can choose a c in (fTB")) (see i j4.4l below) such that the following 
statement holds. 

Proposition 17 (Analyticity in X). XB^-p,i,mX-B is analytic with respect to X on 

the compact set Bf x D. E3 

The proof is given in ">\4:A\ 

As a corollary, we have the following regularity property of £/ m and 1i, m Y . Let 

(26) S- e = {p:lmpe[-e,uj-e), Rep > 0} 



Corollary 18. Forpi € §_ e , define X\ = (p]/ 2 , Zij, where Z\ = cxp [inb/ (2\/—ip\)\ . 
Then, r £i^ m Yo and £; iI7l are analytic in X\ on the compact set Bj x D. 

Proof. Note first that Propositions [8] and [TT] and the relative arbitrariness in the 
choice of p c and e imply that Xb^^Xb is analytic in p in a neighborhood of p = inu), 
n£2\ {0}. Since for large |n|, Xb^/sXb = Xb^oX-b, its expression as an integral 
operator involving G in (|4T?|) (see Note (|28|) as well) implies a lower bound of the 
analyticity radius independent of n. For sufficiently small e for any n € Z, including 
n = 0, then, analyticity of Xb^^Xb in the expanded variable 

Z7T& 



inui, exp 



follows in the domain 



\P- 



into 



|l/2 < 



exp 



(2y — i(p — inui)) 



iixb 

(2y/— i(p — inu))) 



< 1 



^ As usual, by analyticity in a compact set, we mean analyticity in some open set containing 
the compact set. Analyticity in xD of course, implies that XB^-/3,l,mXB is given by a convergent 
double series in p 1 / 2 and e"5X" . 
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(since Proposition [T71 gives analyticity X 6 ED^ x D.) Analyticity of (£; iTn in Xi G 
xD now follows since £i, m is the norm limit of analytic operators (the restrictions 
of <ti : m to the subspaces of H with zero components for \n\ > N). (See Proposition 
[Tl]for the necessary estimates of decay in N.) The analyticity of 1i t rnXo follows 
from its definition. | 

3.6. Compactness. 

Proposition 19. <£i, m is compact in H (cf. Note\^) for p\ G H. 
The proof is given in i j4.5l 



3.7. Step 3. The Fredholm alternative. We can now formulate the ionization 
condition using the Fredholm alternative. 



Proposition 20. // \2J$ has no nontrivial solution in H for p\ € §- e , then 
(I — £i tTn )~ exists and the system ionizes (cf. 

The first part is simply the Fredholm alternative. Ionization follows from the 
following proposition. We recall that y(pi + inu, x) = y n are the components of Y. 



Proposition 21. Assume \24\) has no nontrivial solution when pi G H. Then, for 
^0 € Pz, TO (C^°(R 3 )), the solution Y £% to Wfy is analytic in pi G H\ {iuj'L} and 

analytic with respect to X\ in Dt x D. In particular Y is bounded at p\ = 0. These 
properties imply sufficient regularity and decay of y(p, x) so that the integration 
contour in L~ x y can be taken to be iR. By the Riemann-Lebesgue lemma, P(t, B) — > 
as t y oo . 

The proof is given in ij4.6l 

3.8. End of proof of Theorem [2j It only remains to make the connection with 
Floquet theory. This is done in 



3.9. Proof of Theorem \S[ Eq. ([M)) . restricted to B, follows from the homoge- 
neous system w — £w. Multiplying (|14[) by W^, summing over n, and integrating 
over B„, where a G (l,a], we are lead to a nonnegative definite quantity involving 
w n \dBa being zero for n < 0. Details are given in ? 14.71 

3.10. Proof of ionization for spherically symmetric SI, Theorem [TJ We 

consider the case V(t, x) — 2f2(r) sinwi, corresponding to Sli = — iVL and f2_i = iQ, 
(il is real valued). The proof in the slightly more general case 2f2(r) sin (u>t — 9) 
amounts to replacing t by t — 9/lo and tpo(x) by ip(Q/uj,x) in our proof. Recall 

;i l We obtain by projection of ([23]) to P; iTO (i 2 (R 



(27) Y = Y Q + € lm Y 
The homogeneous equation associated to ([27]) is 

(28) to = € lm w, wen 

The Fredholm alternative applies and (|2"?| has a unique solution in T~L iff (|2"5|) implies 
w = 0. 



(3)As discussed, it suffices to show ionization on a dense subset of initial conditions. 
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Note 22. By separation of variables in spherical coordinates, we see that £i, m can 
be defined in the same way as £, replacing Ap in (IT71) by 

d 2 2 d 1(1 + 1) b 

(29) ^ =-—--— + v ; ^p 1 +nw-i/3XB 

ar^ r ar r z r 

and the associated differential-difference systems are obtained by replacing — A — 

, . . , d 2 2 d 1(1 + 1) b 
b/r W ith - — - -— + - -. 

ar z r dr r z r 

Clearly if there exists a nontrivial solution ro£?iof ([25)1 , then, again by elliptic 
regularity (see Proposition [SJ, v defined by w — yi, m v(r), where 3^,m are the 
spherical harmonics, is a nontrivial solution to 
(30) 

■Ap tr v n = -id {v n+1 - u n _i) - if3X B v n implying Ao, r v n = -iO, (v n+1 - u n _i) 

Proposition 23. If v satisfies H30\) , then there exists n > such that either (i) 
v n (l) 0; or (ii) v n (l) = 0, 6ui f^(l) 7^ 0; Zei no 6e the smallest such n. By 
homogeneity, we can assume that v„ (l) — 1 in case (i) and v' no (l) — in 
case (ii) (we use the positivity of VI). 

The proof is given in H4.10I 

Definition: We define r to be or 1 in case (i) and 1 in case (ii) respectively. 

3.11. Asymptotic behavior of v n in (1301) as n —> —00. In view of Proposition!!?!! 
we see that (f3U)) holds the necessary ionization information. 

3.11.1. Notation. Let 

(31) S(r):= / VW)d P (re (0,1)) 
By assumption SI > is smooth and then so is s. Let 

(32) n -k, fi = O(0), fi' Q = S1'(0), S = s(0), a = , C = a/cr 
Denote 

^0(0 := ^1 e<C 1/2 ^+i/ 2 (C); Go(C) = ^C 1/2 / i+ i/ 2 (C) 

where i^+1/2 and //+1/2 are the modified Bessel functions of order I + 1/2. It 
follows that for small (, 

H (C) ~ 2- l C l {2l)\/l\ 
Let H (Q; fc, I) be the unique solution of the integral equation 

(33) H((;k,l)=G (0 [ ' e~ 2s G Q {s)K(H + AT 1 H)(s)ds 

Jo 

K 



flo(C) / e-^Hois^Ho + k^H^ds 

Jka 



for £ e [0, ka], where the operator TZ is defined by 



14 



O. COSTIN, J. L. LEBOWITZ AND S. TANVEER 



Define 

H(Q = H((; k, I) := H (() + k^H^; k, I) 
It can be checked that H satisfies 

(34) H - 2 { 1_ 2fcc? + 4^ + 2fcJ H + v ^^~<fcJ 
with the following asymptotic conditional 

(35) (0 ^ 1 + + ^a l0gC + ° (rkT' (C ^ ^ °°' C ^ ka) 

Remark 24. fij F(C; k, I) ~ y^e c C 1/2 ^+i/2(C)(l + o(l)) as fc -> oo 

(ii) From the expression \3S\) for H it is seen that as £ — > we have H(C', k, I) ~ 
const. for I ^ 1 and H ~ const. + const Cl°gC / or I — 1- ^ 7 ' or t — or 1, H 
is less singular than H at £ = 0. 



Define 

a; /" r s(s)ds 



36 m fc (r = Y— ^ — exp 

(2k + r)\9.l(r)H(ak) 



4Jl y/Q(s) 



S 2k+T 



(2k 



:F k {r) 



Note 25. From standard properties of the modified Bessel function i^+1/2, it 
follows that for large enough k, H(akr) is continuous and nonzero for r G (0, 1] 
and that as r — > 0, H(akr) is singular as r~ l . Therefore for any k sufficiently large 
it/. = r m& has a finite limit nonzero limit as r — > + . 

Definition 26. With no as in Proposition \23l we define hk(f) by 

i k 

(37) w„ _ fc = —m k (r)h k (r)y^ m 

r 

Theorem 4 (Behavior as k — > +00 (i.e. n — > —00) ). For any sufficiently large 
k, Ufe := r l mk(r) is continuous in r G (0, 1] and Uk(r) — > const ^ as r — > + . 
Furthermore, if there is a nontrivial solution to i30\) . then there exists a subsequence 
kj — > 00 such that for any r G [0, 1], 

(38) lim h k ,j(r) = 1 

3 -too 

The first part follows simply from Note [25l The rest of the proof is given in §4.111 
Proposition 27. There is no nonzero solution of \28\) in %. 

Indeed, Theorem 0] shows that otherwise (r l+1 v n )(0) = m n 7^ for a subsequence of 
n < 0. This implies that the corresponding w n (x) ~ Tn n r ~ 1 y^ m for r = \x\ —> 0. 
This singularity is incompatible with w n G H 2 , (see Proposition [5]) . Thus there is 
no admissible solution of the homogeneous system and the first part of Theorem [l] 
follows from Theorem [2] (i) . See also the remarks in ^5. 101 The result on the decay 
rate follows from the type of essential singularities of y for p G iutZ; see §4.61 



< 4 >As is common, the notation O (a, b) = O (\a\ + \b\); similarly O (a, b,c) = O (\a\ + \b\ + |c|) 
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4. Proofs of intermediate steps 

4.1. Proof of Proposition |8j As mentioned in Ap and At are adjoints of 
each-other. They are furthermore densely defined and hence (see, e.g Theo- 
rems 5.28 and 5.29, p. 168), closed. Once we show that Ap(p) is invertible in T>, 
analyticity of Tip in O \ (£ Pc U£_ e ) follows (the spectrum of the closed operator 
H c — iftX B {r) is a closed set). Analyticity holds wherever Ap is analytic, [3T], vol. 
1, Theorem VIII.2, p. 254). 

(1) Eigenvalues. We first show below that no ip S iD is an eigenvalue of He — 
i(3X B (r). Assume we had Apip = 0. If /3 — 0, Aoip = implies ip £ <r p (Hc), but, 
by construction, these values of p correspond to the region where (3 ^ 0. So we can 
assume f3 > 0. Then 

(39) (ip, (-A - ip - br- 1 )^) + (ip, -i(3X B ip) = 
Taking the imaginary part of (|39p we get 

(40) Rep (ip, ip) + c(ip, X B ip) = Rep (ip, ip) + c(X B ip, X B ip) = 

If Rep > this immediately implies -0 = 0. If Rep = we get X B ip = 0. But 
X B ip = implies = Apip — Aoip. In spherical coordinates the equation Aoip = 
becomes a system of ordinary differential equations 

( d2 2 d Kl + 1) , i A , 

(41) {-^-rdr + ^ 1 - br - 1 -ip)^. m = 

Since X B ip = 0, the solution of (|4"Tj) vanishes identically on [0, a]; but then, by 
standard arguments the solution is identically zero. 

If Imp ^ [— e,p c ], with p £ T>, then _4,g = Ao, and we are, by construction, 
outside the spectrum of Ao, and thus Aoip = implies ip — 0. 

(2) TTie range of Ap is dense. Indeed, the opposite would imply W\ Ker (A % ) ^= 0, 
which leads to the same contradiction as in Step 1 (note that A*p is simply Ap with 
the signs of j3 and Rep changed at the same time). 

(3) For any p G T> there is an e > such that \\Aptp\\ > e||'0||. 

(a) If Rep > and \\ip\\ = 1, then 

(42) \\Apip\\ > \{Apip,ip)\ = \(A ip,ip) -ic(X B ip,X B ip)\ > \Rep{i>,ip)\ > Rep 

(b) Let now Rep — 0, and assume Imp is between two eigenvalues of —He, the 
distance to the nearest being 5 > 0. To get a contradiction, assume that \\ipj\\ = 1 
and H-A^VjII = £j -> 0. Then 

(43) 

£j = \\ApiP 3 \\ > \{Ap^ 3 ,iP 3 )\ = | (Aoipj , ipj) — ic(X B ipj , X B ipj) | ^ c|(X b ^>Xb^)I ->■ 
thus X B ipj — > 0, and by the definition of Ap and Ao we get 
(44) HA^II^O 

which is impossible, since our assumption and (|44p imply noninvcrtibility oiHc — ip 
while ip is outside the spectrum of Hq- 

(c) In the last case, Rep = 0,lmp G a p (—Hc); then if we assume there is a 
sequence ipj, \\ipj\\ = 1 such that ||»4a?/>j || — > as j — > oo we get 

(45) 



( 5 )[28], p. 267 



16 



O. COSTIN, J. L. LEBOWITZ AND S. TANVEER 



Since H-Ao^jH < 11-4/3^11 + c\\X B ipj\\, (HSJ) implies HAV'jII -> 0. On the other 
hand, with P the orthogonal projection on the finite dimensional eigenspace of Hq 
corresponding to the eigenvalue ip, we have AqP = =>■ Ao = Ao(I — P) and then 
since Aoipj —t 0, 



But by definition Ao is invertible on (7 — P)L 2 (U. 3 ) and (l4"6"| then implies ||(J — 
P)ipj\\ -> 0, i.e. Pipj - tpj -> 0. Since \\ipj\\ = 1, \\Pipj\\ -> 1. Then P^j is a 
bounded sequence in the finite dimensional space PL 2 (R 3 ), hence we can extract 
a convergent subsequence, which we may without loss of generality assume to be 
Ptpj itself, Pipj — > ip, HV'II = 1, and also ipj — > Pipj —> ip, thus Pip = ip. Therefore, 
Aotp — AoPip — 0. Also, since multiplication by cX B is a bounded operator we have 
cX B ip 3 -> cX B ip = 0, since cX B ipj -> 0. Therefore, \\Apip\\ < \\A tp\\ + \\cX B ip\\ = 
in contradiction to the absence of eigenvalues. 

(4) Definition of the inverse. This is standard: we let ip € D(Ap), Apip = <j> 
and define %\p<p = ip- This is well defined since Apipi — Apip% entails, by Step 1, 
ipx = ip2- By Step 2, is defined on a dense set. By Step 3, for any p there is 
an e > such that \\9\p\\ < e _1 . Thus d\p extends by density to L 2 (R 3 ) and by 
construction Affipcp — (p whenever *Rp<p e D(Ap). Conversely, if (p G -D(^4.,s), and 
^0 = u then SH^u = <p entailing D\pAp(p — (p on the dense set D(^). 

For the regularity of in i, we first note that if we define Q = (I — A)" 1 , we 
have the following identity: 



It is clear that if <p € L 2 (R 3 ), Q<P € H 2 (R 3 ) and so (b/r ~ i/3X B + ip + I) Q<p € L 2 . 
Therefore, from fl47|, % : L 2 (E 3 ) -> ff 2 (R 3 ). 

4.2. Proof of Proposition 1151 The shift operator 5, defined by (5y)j = 2/j+i, is 
quite straightforwardly shown to be bounded in T-L: the proof of Lemma 27 in [TS] 
goes through without changes. By the second resolvent identity we have 



Since —A — br 1 is self-adjoint, we have by the spectral theorem, for some C > 
independent of p, 



bounded (cf. Lemma fT2"j) . we have \\£\\ n C 2 (l + (Repil) -1 . Thus \\<t\\ n is small 
for large Repi and therefore (I — €)Y = has unique solution Y € % and the 
proof follows. 

4.3. Proof of Proposition ITT! 

Proof. The estimate ||Xb9^oXb|| = 0(p^ 1 ^ 2 ) is shown right after the statement of 
Proposition [TT1 

We now consider the analyticity of 9\p in an open set on the imaginary p axis 
for Im p < — e. There, = and Xb^-^Xb = Xb^oXb is manifestly analytic from 
its representation as an integral operator, whose kernel G is given below (see [21)] 
and Appendix ^5.41 for details). 



(46) 



\\A (I - P)ipj\\ ^0 



(47) 
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With k = ^/vp (using the principal branch of the square root), and v = b/ (2k), 

(49) G(x,x';k) 

_ ik(r) - £)I(-ik£)J(-ikr)) - k 2 £ri[I(-ikt)j(-ikri) - J(-ikrf)i(-ikCj\ e^tf-H) 
~ r(l - iv)T(l + w) 4n\x-x'\ 

where 

(50) £ = \x\ + \x'\ + \x-x'\, r}= \x\ + \x'\ -\x- x'\ , 

pioc pioo 

(51) I{zi)= e- Zlt r w {l+t) iv dt, i( Zl ) = - e-^H^il+tfdt 

Jo JO 

J(z 2 )= f e z *H- w {l-t) w dt\ J(z 2 )= f e z ^t x - iv (\-t) iv dt 
Jo Jo 



Further properties of function G are discussed in ^5.41 

Note 28. Note that (03) still holds for p £ iR+, with k = y/ip, with the choice 
arg/c — tt/2 for p e iR + , and with the upper limits ioo in (|5ip replaced by +oo. 

4.4. Proof of Proposition I17[ The function / = ^/3,i,m^B9 is the solution of the 
equation 

d 2 2d 1(1 + 1) b n2 . \, 

1 \ 2 -ic )f = g; r^a 



dr 2 r dr r 2 r 
( d 2 2d 1(1 + 1) b x2 \. n 

-TT ~-^-+ 2 -~ + A 2 / = 0; r>a 
\ dr z r dr r z r / 

(52) 

such that / decays at infinity, is regular at the origin and C 1 at r = a. We note 
A = \/—ip is in the closure of the fourth quadrant for Rep ^ 0. We let a = \/X 2 — ic, 
Ki = b/(2a), k = b/(2X), p — 21 + 1 and define (in terms of the Whittaker functions 
m and 

(53) 

mi(s) := s~ 1 M Kl ^/ 2 (2as); ioi(s) := s~ 1 W Kl ^/ 2 (2as); w 2 (s) := s^ l W K ^ /2 (2\s) 

For r > a we have / = Hw 2 (r) since r^ 1 dJl K fl ^ 2 (2Xr) grows with r as r — > oo. For 
r^awe must have 

(54) / = Arm + jo 

where, using standard results about the Wronskian of 9Jt and 2U, see [9], pp. 25 
and pQ, pp 505, 508, we have 



2ar (l + p) 

/ X[o,o](s)s 2 TOi(s)5(s)ds + mi(r) / s 2 ioi(s)X[ ,a](s)fl , (s)rfs 



( 55 ) r (i , i 







( 6 )See [9], pp. 60 eq. (1) and pp. 63 eq. (5). 
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The integral representations of the functions DJl and 2U < - 6 ' entail immediately that 
the functions /o, / and 9Jt rei ljU 12 (2ar) depend analytically on A for small A. Conti- 
nuity of / and /' at a ^ 1 imply that A defined in ([51)1 . is given by 

^ ^_ /oO)?4(a) - f Q (a)w 2 {a) 



m[(a)w2(a) — mifajw^a) 

In §5.51 it is shown that that A is analytic in (A,exp [i7r6/(2A)]) in a domain cor- 
responding to A small in the closure of the fourth quadrant, if a and c are cho- 
sen large enough. It follows that resolvent 9\p, m ,n is analytic in X for X = 
(y/p,cxp [iirb/ (2y^ip)]) e Bf x B for small e. 

4.5. Proof of Proposition 1191 By adding and subtracting 1 from Ap and using 
the second resolvent formula, whenever everything is well defined, we have 

(57) XbA^Xb =: XbK/sXb = X B (-A + 1)- X X B 

- XzdXpi-br- 1 - ipX B - 1 - ip)(-A + 1)-% 
The Green's function for — A + 1 is 

(58) 0(2,9)= 1 ■ e-'*-'" 

An\x — y\ 

Now if ||</>j||L 2 (B) ^ 1 then the functions fj — (— A + 1) _1 Xb0j are seen by straight- 
forward calculation to be equicontinuous on the one point compactification of R 3 . 
A subsequence, without loss of generality assumed to be the fj's themselves, con- 
verges in L 2 (R 3 ) as well (to a function with exponential decay, since there is a 
Si > small enough and independent of j so that e Sl ' x '(—A + 1) -1 Xb^/) is also 
equicontinuous on the compactification of R 3 ). In particular, Xb(— A + 1) _1 Xb is 
compact. 

Now fj converge in the sup norm with weight e Sl ^ , and thus (— br^ 1 — iflXs — 1 — 
ip)fj converge in L 2 (R 3 ). Since dKp is bounded, compactness of Xb9*sXb follows. 

By Proposition [TlJ and the previous argument, £ is a norm limit of compact 
operators (the truncations of £ to the subspaces of % with vanishing components 
for |n| > N). Therefore, <£i : m = Pi, m £ is also compact. 

4.6. Proof of Proposition 1211 and final estimates for Theorem [1} If 

has no nontrivial solution for any p\ € H, then compactness of £/ !m implies that 
(I — £/ !m ) _1 exists. Lemma [THl Corollary IT51 and Proposition ITT1 give the analytic 
and continuity properties of £; jm and 2^ m Yo- Analyticity of (J — £/, m ) 1 in X\ for 
X x e B+ x B, follows in a standard way from analyticity of €-i, m and the second 
resolvent formula, 

(59) A- 1 - B- 1 = B- l {B - A)A^ 

(see £I5.3[) . The same resolvent identity can be applied to show analyticity of (I — 
£l,m) with respect to p\ in a neighborhood of 

H \ {(V + U {t-e + i^L) U (X e + iwZ)} 

Hence, the solution F = (I — £; jTO )~ is analytic for p\ G H \ {iujZ}, since 
e, /3 and p c are artificially introduced parameters the value of which cannot affect 
Y, since Y is independent of these choices (see Remark UM ) 
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The function y(p;x) — y n (p%,x), with p — inuj + pi, is analytic in p for p £ 
iM. \ iuiZ and by analyticity of Y in Xi, boundedness at p = iu)7L follows. In 
particular, as p — > into from the right half-plane, y(p, x) is analytic in the extended 
variable 



(60) 



(p — inu;) 1 / 2 , exp 



2yJ—i(p — into) 



The regularity properties of Y in p and the decay properties in \n\ of its com- 
ponents y n for large |n|, simply stemming from Y £ Ji, imply that y(t,x) can be 
expressed as an inverse Laplace transform of y(p, x) on iR. We now show that 



P(t,B) = WM^e-'+y^WL^B) < ^- 2t \\^o\\h {B) +2\\y(x,t)\ 



L*(B) 



as t — v oo 



We note that 
(61) / dx\y(t 7 x)\ 2 



dx 



OO /-OO 



e »t(s s )y(is,x)y(is' ,x)dsds' 



OO — OO 



y(is + is' , x)y(is', x)dx 



ds' \ ds 



So, in order to show ionization, it suffices from Ricmann-Lebesgue Lemma to show 
that 



y(is' ', x)y(is + is, x)dx 

is in L (ds). This follows from Cauchy-Schwarz, since 
(62) 

/•oo f poo r /• 

(is' ,x)\\y(is' + is,x)\dx ds' )■ ds < 



ds' 



oo w — oo 



\y(is'i -)\\L^(B)ds' 



However, 



(63) / \\y(is' ,-)\\ L * m ds' = 



\\yn(iq,-)\\L*(B) 



dq 



<C V(l + \n\)^ 3 \\y n (iq, OllkfB)*? < C / \\Y(iq, -)\\ 2 H dq < oo 

since Y is bounded in p\ = iq for q £ [0, w]. 

Since y(p\ + inuj,x), n £ Z, is analytic in the variable (|60|) . standard stationary 
phase analysis (see Appendix ^5.71 ) shows that y(t,x) = 0(i~ 5 / 6 ), and hence 
P(t, B) = 0(£- 5 / 3 ) ast-^oo. 



4.7. Proof of Theorem |3j Since (fHij) (restricted to B) follows from the homoge- 
neous system w = €.w (see also Proposition [8] for the necessary regularity), we look 
for a nontrivial solution of ([T4")) in H. We multiply (fT4"l) by u; n , integrate over the 
ball B a (of radius a £ (1, o]), sum over n (this is legitimate since w £ H) and take 
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the imaginary part of the resulting expression. Noting that 



(64) ^2 ttj(x)w n -jw n = n_ 

j',neZ j,ne'. 



jW n -jW n — * L f w n+j w n 



so the sum (|64p is real, we get from (fbf)) 



(65) = Im 



E / \w n (x)\ 2 dx + / dxw n Aw n J 
Rep!^2 \w n (x)\ 2 dx+^ l^2w n V 



— w n Vw n \ ndS 



It is convenient to decompose w n using spherical harmonics. We write 

W n = Rn,l,rn(r)yl n (0,^). 
l^0,\m\^l 

The last integral in ([US")) , including the prefactor, then equals 



(66) 



(67) 



■ a 



•} 



a 2 ££W[i? 



where W[/, <?] is the Wronskian of / and g. On the other hand, we have outside of 



Aw„ + br 1 w„ + (ipi — nuj)w n = 
and then by (|66[) . the R n j, m satisfy for r > a the equation 



(69) 



n 2 . , Z(Z + 1) 



where w6 have suppress6ci th.6 subscripts. Let (jn,i.m = rRn,i, m - Then for the g n ^. 
we get 

' 1(1 + 1) 



(70) 
Thus 
(71) 
and 
(72) 



.9 



— ipi + nui — br 1 



.9 = 



r 2 W[i?,i?] = W[g,g] =:W n 



Multiplying (I7n)) by g, and the conjugate of ([70)1 by g and subtracting, we get for 
r > a, 

(73) V^ = -i(pi+pr)|ff| a = -2i| ff | 2 Hep 1 
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Remark 29. Direct estimates using the Green's function representation \49( j imply 
that 

(74) w n {x) = \ + _^_ (c n {0, 4>) + 0(r -1 )) asr^oo 

with c n (9,4>) independent of r and with 

(75) ^ 

n n = \j — ipi + ulu (when Repi > ; K n is in the fourth quadrant when n < 0) 

(i) We first take Repi > 0, to illustrate the argument. Using ([74]) we get 

(76) ff ~Ce- K " r r~5^(l + o(l)) as r -> oo 

There is a one-parameter family of solutions of (170)) satisfying (1751 and the 
asymptotic expansion can be differentiated [42] . We assume, to get a contradiction, 
that there exist n < for which g = g n ^ 0. For these n we have, using flT6"j) . 
differentiability of this asymptotic expansion and the definition of n n that 

(77) i lim \g n \- 2 W n = -Im K „ > 0. 

It follows from ([75)) and (177)) that W n /(2i) is strictly positive for all r > a (by 
monotonicity and positivity at infinity) and all n for which g n =/= 0. This implies 
that the last term in (|65p is a sum of nonnegative terms which shows that (|65p 
cannot be satisfied nontrivially. 

(ii): Kepi = 0. For n < 0, we use Remark (and differentiability of the 
asymptotic expansion as in Case (i)) to calculate W n in the limit r — > oo: W n = 
2i|c„| 2 |K„|(l + o(l)). Since for Repi = 0, W rl is constant, cf. ((75)) . it follows 
that W„ = SilcnplKnll^nl 2 exactly. Thus, (|65)) cannot be non-trivially satisfied, 
implying that 

(78) w n (x) = for all n < and |x| = r = a € (1, a] 

For a > r > 1 (where V(t, a;) = 0) we have Ow n = 0, where D is the elliptic 
operator —A — b/r — ip\ + nuj. The proof that w n (x) = for r > 1 then follows 
immediately from ([78")l. by standard unique continuation results j23j[17j (in fact, D 
is analytic hypo-elliptic). Sec also Note [5] 

4.8. Connection with the Floquet operator. It is easy to check that the dis- 
crete time-Fourier transform of the eigenvalue equation for the Floquet operator, 
Eq. ([5]), Kv = (j>v, with p\ = i<fi, coincides with (|14|). the differential version of the 
homogeneous equation associated to (J23J. Now, ((78)) shows that a solution of ((14)) 
is an eigenvector of K. 

In the opposite direction the existence of a Floquet eigenfunction entails failure 
of ionization since it implies the existence of a solution of ^ for which the absolute 
value is time-periodic. 

4.9. Differential equation for w. We seek to show that the only solution to the 
homogeneous system 

(79) w = £i m w 

in the space W is w = 0. Since w is piecewise C 2 (see Note 03), ((79")) implies that 
the components of w = {3^,ro r " _1 <7n(f)} neZ satisfy the differential-difference system 
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(see Note El): 

(80) ~^2 9n ~ (^ br ^ 1 ipi + l(yl ^ 2 — \ g n = id (g n+1 - g n -x) 

First, we notice that for n < 0, Theorem [2] implies that g n (r) = for r ^ 1. Thus 
g n (l) = 0, g' n {l) = for all n < 0. 

4.10. Proof of Proposition I23[ The gist of the proof is that contractive mapping 
arguments show that if the statement was false then the solution would vanish. 

Lemma 30. IfY^O, then there exists some n a ^ so that either 9n (l) =/= 
or g' n (1) =/= 0. (As before, in the sequel, we shall define no to be the smallest such 
integer.) 

Proof. To get a contradiction, assume the statement is false. Since the functions w n 
are in the domain of A (see Note [5]), then, in particular, for any n, g n is continuous 
in r. Thus, the set Z n := {r : g n {r) — 0} is closed and so is the (possibly empty) 
left connected component of 1 in Z ril call it K n . Let 

k = p| k u 

Assume to get a contradiction that K is nonempty: let then K = [a, 1]. If a = 0, 
then Y = since g n (l) = 0, ^(1) = imply g n (r) = for r > 1. Then Y ^ 
implies a > 0. We first take < a < 1. We write the differential equation for g n (r) 
in integral form and use the conditions g n {a) = = g' n {a), since g n vanishes on 
[o,l]: 



(81) e^ r g n (r) 




I _ g— 2\/nu{s— r)l 



W + 1) • 

- i Pl + V(s) 



s* 

Consider the Banach space of sequences 
in the norm 



j n (s) - ifl(s) (g n -i(s) - g n+ i(s)) } ds 



sup 

n£Z,r^ \a—e,a~\ 



lulr 9n(r) 

It is easy to see that the rhs of ([ST]) is a contractive mapping if e is small enough 
and then g n (r) — for r E [a — e,a] contradicting the definition of a. The same is 
true if a = 1, since g n (l) = and g' n (l) = would imply, with the same proof as 
before, that g n = for r e [1 — e, 1], for some e > 0, contradicting the definition of 
a. I 



4.11. Proof of Theorem |4| For a heuristic discussion see 

The proof is by rigorous WKB. The fact that there are two competing potentially 
large variables, k and 1/r makes it necessary to rigorously match two regimes. 
First, note that (f37|) implies 

(82) 9n -k(r) = i k m k {r)h k (r) 

We need a few more preliminary results. 
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Lemma 31. For any e\ > 0, there exists C3 > independent of k and t\ so that 
for k ^ fc = C'3 e i 1 > an d f or r € [ei, 1], 

fk \ 1/2 

(83) sup \h' k \^C 4 k -M 

ei <r<l V fc / 

where C4 is independent of ei and k. 

The proof of Lemma [3T] is given in i)4.13l 

Definition 32. For fixed e, we define L t — aC^ (2C4C3/ 'e) 2 , withC^ andC^ defined 
in Lemma lSH and C = afcr, where a is given in iSty) . We will take e small enough 
so that L e > C3CK. 

Finally, in what follows, c + is a positive "generic" constant, the value of which 
is immaterial. 

Lemma 33. For e > small enough and kar — £ £ [Z e , fca], we have 

(84) |/i fe (r) - 1|< e 
The proof of Lemma [33] is given in ^4. 141 

Definition 34. Let h k (() = h k (C/(ak)). 

Lemma 35. For any small e > 0, there exists a subsequence S = {hkj }jeN that 
converges to a continuous function h for £ £ [0, L e ] . For the limiting function h(£), 
we have \h(() - 1| < 4e for ( £ [0,L e ]. 

The proof of this proposition is given in i)4.15l 

Proposition 36. For any r £ [0, 1], lim^oo hk,j(r) = 1. 

Proof. From Lemma [33] and Lemma [3jJ] it follows that for any r £ [0, 1] and any 
e > we have lim^oo \hf tj (r) — 1| < 4e. | 

The proof of Theorem 0] now follows from the definition of hk in (|36|) . Remark [24| 
Note [5] and Proposition l36l 

4.12. Further results on g n<) -k and /ifc. 

Lemma 37. for any j, k £ N U {0} we have, at r = 1, i.e. at s = 0. 

Proof. In case (i) (corresponding to r = 0), note that ([50)1 may be rewritten, cf. 
ID, as 



(85) (9no-k)s3 — 2Q3/2 (9n -k)s + ~^~9no-k — i (ffn -fc+l — 9n -k-l) j 

where 

(86) Q fe = - + (fc - n )w + ip x - M+il 

r ?' z 

Since <7„ _fc(l) = = fc (l) for all k ^ 1, while <?n (l) = 1, the statement follows 
from (|85[) for any ^ j ' ^ 2, if 2k j. Assuming the statement holds for some 
j ^ 2 for 2fc ^ j, we prove it for (j + 1) for 2fc > (j + 1). 
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2:21; 



Taking (j — 1) derivatives in s of (f85f) at S = 0, we obtain 

d j+1 9n -k _ . &- x . 

foj+i ~ l d^f=l 9n o-{k-X) ~ *^n<?„ -( fc +i) + L 

where L is a linear combination of derivatives of g no -k up to order j, which are all 
zero since 2k ^ (j + 1) > j. The first two terms on the rhs give a contribution of 
M fc <5(j-i), 2 (fc-i)+0 since 2k ^ + implies 2(fc- 1) ^ (j - 1) and 2{k + 1) > (j-l) 
completing the inductive step. 

In case (ii) (corresponding to r = 1): since <?n (l) = and g no - k (i) = = 
Sno-fcW ^ or a ^ k ^ 1, it follows from that g^ _ k = for all k ^ 1 implying 
the conclusion for j = and j = 1. By taking an additional derivative of (|55|) with 
respect to s and evaluating at s = 0, we obtain 

d 3 g no -k ., 9 c4 (l) 

— 102.2k^-g no s=0 — *<>2.2fc , . r 

9s 3 9s ~\A2(1) 

so the statement holds for j = 2 and any fc with 2fc ^ j. The rest of the proof is 
very similar to that for r — 0. 1 

Let ip2,k be two independent solutions of 

(87) £ fe V = ; and W k = ^i, fe (r)^ fe (r) - ^(^(r) 
where 

(88) = ^" + g fc ^ 

From the form of the equation we see that W k is independent of r. 
Lemma 38. For n = hq — k, k ^ 1, i/ie system I180\) is equivalent to 

(89) g no - k (r)=i Q(s) (g no - k +i(s) - g no _ fc _i(s)) G k (r, s)ds k^l 

J r 

where 

(90) G k (r,s) = Wfc- 1 [V>i,*(r)V2,*(«) - V'2,fc(r)V'i,fc(s)] 

Proof. The proof simply follows from variation of parameters, the two boundary 
conditions at r = 1 and g„ _ fc (l) = <^ _ fc (l) = 0. | 

Definition 39. Define 

5 

(91) j k = [C k m k - VLm k _i] 

m k 

Lemma 40. For k > 1, i/iere exist constants C\, Ci and c*, independent of k 
so that for any r 6 (0,1] rac Ziaue < c* . for r > r, we ftawe < 
d/(fcr 2 )+C 2 

Proof. In the Appendix, (|253p . we obtain an explicit expression for jV Routine 
asymptotics for large k in different regimes of r £ (0, 1], discussed in the Appendix 
M5.81 show that k 2 j^ + kjjp = 0(1) in all cases and hence j k = 0(1). In fact, 
as r — > and fc — > oo with £ = /car = 0(1) fixed, we have j k — >■ 5(C), where 
g(() is bounded. Also taking the r- derivative of j k for r — 0(1) not small, we 
get jjL(r) = 0(1). When r <C 1, the asymptotics in the regime | <C r <C 1 gives 
j k = O (C _1 ) = O (l/(fer)). Since the asymptotics is differentiable, we have j' k (r) — 
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O (l/(kr 2 )). Finally, we look at C - 0(1), C > 1- Since £j k = k^ Jk ~ kg'(() 
where C 2 .°'(C) is bounded for all (, it follows that \j' k (r)\ < C\j (fcr 2 ) + C 2 for 
r>l/k. | 

Lemma 41. For k ^ 1, h k (r) defined in ilffiij) satisfies the system of differential 
equations: 

(92) /i fe + 2/i fe — * + 1 /ifc = S2 hk-i{r) H /ife+i(r) 

m k \ m k 5 J V m fc m fe 

and i/ie system of integral equations &89\) is equivalent, for k ^ 1, to 

(93) fe fc (r) = " (8)m ;- l(8) G t (r, 8 )> it - 1 ( 3 )<fa 

1 ft(s)m fc+ i(s) 

^ Gk{r, s)h k+1 (s)ds := A k h k -i +rt k h k+ i 

m k (r) 

Proof. This simply follows by substituting g no -fc( r ) = i k nT-k(^)h k (r) into (|80jl and 
(|89|) . and using 

ml' £/tm fc fim fe _i j fe 

\- Wk = = 1 , 

m k m k m k 5 

in turn a consequence of Lemma 1401 | 

Remark 42. Let now r £ [e, 1], where e ^ C^A; -1 for sufficiently large C2 indepen- 
dent of k. It is convenient to rewrite A k and T-L k in (pl3l in terms of s (see (|3. 11.1(1 ). 
Furthermore, changing the variable of integration from s to t — s(s)/s(r), we obtain 

(94) [A k h k - 1 ]{a) = (2k + T)(2k + r-l) [ t 2k - 2+T T k (s,t)h k ^( 5 t)dt 

Jo 

where, using (f3T>)) . we get 

(95) T k {s,t) = G k {r{s),r{st)) 

sF k {r{s)) 

and 

(96) [H k h k+1 ](s) 5 



(2fc + 2 + r)(2fc + l + T) 

x / vW^ f w?! ) G i (r(i),r(tt))/ lH i(it)iit 

In evaluating ,4fc for large it is useful to calculate the Taylor expansion of T k (s, t) 
and its s derivative at t = 1. To do so, we first note that 

(97) ^ = (_^fe^_W Gt(r , r0 _^^ (ri/) 

di V 2il{r')F k {r) F k (r) J F k (r) dr' 

where, to simplify notation, we wrote r(s) = r and r(ts) = r' and used d r 'S(r') — 
-y/n(r'). From and O we get G fe (r,r) = and d r ,G k (r,r') = 1 at r' = r; 
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(|97|) implies 
(98) 



dt 



F k -i(r) 
' F k {r) 



Using (|9T|) . taking an additional derivative with respect to t, using also (j8"6")l and 
to see that d r > r >Gk = —QkGk, we obtain 



(99) 



d 2 T k 



dt 2 



F fc _i(r) / £fi'(r) 



t=i 



F fc (r) I 2fi 3 / 2 (r) ^ v^J^.^r) 



A similar calculation can be carried out for the third derivative. We only write 
down the potentially largest term in the regime kr > C2 (for large k and small r) 
(100) 
d 3 T k 



e 2 ^-iMQ fc M , ( j_\ = eF.Mr) ( koj _ 1(1 + 1 



dt 3 t=i fi(r).F fc (r) ' " ^' fcr " 3 / ^(r)F fc (r) 

Note that if kr is sufficiently large, (j55")) gives 

^ib-i(r) H(a(k~ l)r)H(ak) 



1, 



1 

fcr 3 



(101) 
and 
(102) 



F fc (r) H(akr)H(a(k-l) 



K-i(r) 1(1 + 1) 
F k -i(r) ~ 



2akr 2 



O 



l + 0(k- A r- L ) 



k 2 r 3 



Note also that (|52"j) implies a — 2y/Q(r)/ (s(r)) = O(r) for small r. Including all 
terms that become important when r is small, we note that in the regime when kr 
is sufficiently large, we have 



11.08] T, = (1-/)+(4/i + S] U(!-*) 3 - 



(1-t) 5 



(1-t) 4 (1-t) 3 (l-<) 3 (1-t) 2 (1 -i) 2 (1 



fcr 3 



kr 



k 2 r 3 k 2 r 



(104) 

where 
(105) 



m 

ds 



= [~ fx 



k , M 



1(1 ~ 



4* x ' r 3 J \3 y '' k 

()l (i-t) 4 (1-t) 3 (i-ty a-n- u-/> 2 u- /) 



Q(r(s)) 



(106) 



/2(B) 



Z(i + l)s 2 
417 



(107) 
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When r € [0, e], for e = Cz/k, it is sometimes more convenient to express *4fc in 
terms of £ = fear. For that purpose, we define 



(108) Q(C) = -2fclog 



1 



s(0)-s 



log 



exp 



1 r±j 



dr 



where we recall the relation (f5Tj) between s and r = £/(ka), C S [0, fcae]. A series 
expansion in /c _1 leads to 



(109) Q(C) = C - 



C / w n'(o) \ c 



k \2a 2 4fl{0)c 



4fc 



1 



n'(o) 

afi(0) 



We choose i\ = C^k 1 logfc, for some fc-independent C2 (chosen more precisely 
later). We define Si, dependent or r, so that 

(110) {l-5 1 )s(r)=s(e 1 ) 

From (f3Tj) . it follows that for sufficiently large C2 we have 



(111) 



S t p 1 s ^ > ( 5 + 0iogfc 



s(e) (4fc + 2r)fc ' 
It follows from the definition of Ak in (|93|) that for r <E [0, e], i.e. ( £ [0,kae], 
(112) 

ai\ ^(77(1 - fc" 1 )). 



[A^-i] (C) = 





e -Q(n)+ 










r S (r')" 


2fc- 









+ (2fc + r)(2fc + r-l) 



fi(/)G fe (r,rO Fk ~ l{ ^ h k ^(r')dr 



s 2 ^(r) 



where 77) is defined by 

(113) 0(C,^) = *«G fc (r(O,rfa)) 

while 



(114) 



Oi(»?)C) _ H(ak) 



k 



H(a(k- 1)) 



1 



2k 



1 



- 1\ s 2 (0)n{T 1 /{ka)) 



2k ) s 2 (rj/(ka))n{0) 

' s(rj/(ka)) 
*(C/(*a)) 



while for large /c and 0<£<?7<eia we have 



(115) a 1 ( V ,0=r-l-+ I H 



n 



»n(o) / 



ry+-(C-?/) + 



77^ 77 



Similarly, for fcar = £ G (0, fceia), we define 



(116) 



&ifoQ = g (« fc ) ^^/(fca))^^/^)) / g(7?/(fc a )) \ T 1 



ff(a(fc + l)) s 2 (0)fi(0) 



0(C/(*a)) 
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We then have 
(117) [n k hk+i] (C) 



O(0)s 2 (0) 



a 2 k 2 (2k + 2 + T)(2k + l + T) 



-Q(v)+Q(0 



C 



k) H(Q 



+ (2fc + 2)(2/c + l + 2r) 

VWW)G k (r(s),r(st)) t 2fc+2+T ^ + i(*t)df 

=: [Kgft fc+1 ] + [H^ fe+1 ] 

Lemma 43. for fc ^ 2 and k\ £ {k — 1, fc, fc + 1} we /iave 

(1) If r £ (0, 1) and s G (r, r + 6), where S < min {C^fc -1 logfc, 1 — r}, i/ien 



G k (r, s) 



F k {r) 



5 
9r 



F k {r) 



< c, 



(2) If re (0,1), 5< C 2 fc _1 logfc withr + S < 1, then for s £ (r + 5,1), 



G k (r,s) 



*k(r) 



< 



r 1.//2-1/2 



5_ 

9r 



< 



Li/2+1/2 



Proof. It suffices to find bounds for Gfc(r, s)H (akis) / H (akr) since the other func- 
tions involved are regular everywhere for r,s £ [0,1], see ((55)) . We first consider 
— > +oo. 

It is easily verified that G{(,, r]), defined in (|113p . is the Green's function (see 
) for 



(118) 
and is given by 



T ^ 

-I 

C 2 k 



LO b 

or cm. 



k 2 a- 



[ipi - n uj] 



(119) 



G(C,v) :=kaG k (r(0,r(v)) = 



w 



where "J 2 are two independent solution of = and W — (C)^^ (0 — 
\E f 2(C)^E r i(C) is their constant Wronskian. 

Standard asymptotic results show there exist two independent solutions ^2 
such that for large k, we have uniformly in z £ [0, V ujk] 



(120) 



2'Z! /ttI 



(20! 



^+1/2(2) ; where z 



a 2 k 



c 



ojkr 



(121) 



-L-l 



(2Z + 2)! /ttz 



(1 + 1)! 



The Wronskian is asymptotic, for large fc, to {2l + l)y/uj / y/a 2 k. The expressions 
(I120p and (|121l) may also be used to determine the asymptotics of ^ and 
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Using ([TT9]) . (fT20)) . (fl2Tj) and (H| and the bounds on W, with it = I + 1 it follows 
that 



(122) 



F k (r 



G k (r,s) 



c*\zz I 



'1 1/2 



if ( aJ k z 



[Y h {z)J h {z')-J h {z)Y h (z')] 



where z' — r\sfwj\ 



'uks. A similar bound holds for 
d (F kl (s) 



dr { F fc (r) 



•Gk(r,«) 



We now prove part (1). We break this case up into two subcases: (a) r G 
[fc~ 2 / 3 ,l] and (b) r G [0, k~ 2 ^ 3 ]. In case (a), we note that s G [r, r + <5] implies 
s/r and therefore 1 < z'/z — 0(1). The function H in (|122|) is close to 1 because 
its argument is large. Furthermore, note that y/zY l+1 / 2 (z) and \/ r zJi+i/2( z ) are 
bounded for large z, while they are asymptotic to constant multiples of z~ l and 
z l+1 for small z. Using (|122[) , part 1 of the Lemma follows by inspection in case 
(a). For the case (b), (|122[) further simplifies since 

(123) 



z,z' are small and 



H(k lV /k) 



G*(r(C),rfa)) 



H(kirj/k) 



H(kir]/k) (r] l+1 C 



H(C) y yi " kaH(() kaH{()(2l + l) 

When ( G [log fc, afc 1 / 3 ] and 7/ G [C, C + ctkS], we have 1 < [n/Q L < and therefore 



H(kiT)/k) 



mo 



G k (r((),r(r,)) 



H(k lV /k) 



kaH(() 



The same inequality holds if £ G [0,log/c], since n G [C, (C*2 4- l)logfc] since in this 
regime £~'/ii(£) is bounded and the logarithmic growth in k of terms involving 
rj can be bounded by, say, k 1 ^ 2 , while for small n, n l H(kin/k) is bounded. The 
bounds on derivatives follow in a similar manner using = ka-^. 

Part 2 (which is only relevant for r + 6 ^ 1) follows similarly on careful inspection 
of (|122p . from the asymptotic behavior in different regimes of z and z' . | 

Lemma 44. Let r G (0, e], with £i = \ogk. We choose C 2 large enough so that 

^ = (! - < ™en l^/](r)| < C^'/2-l/2(l - (5l) 2 fe - 2+ r|| / || oo ^ 

c*fc- 3 ||/|U and \£[Alf](r)\ < cfcW/^l _ s^k-2+ry^ ^ c=fc fc" 2 1| / 1| ^ . 



Proof. Consider _4jf, given by (|112|) . We note that s~ 2 fl(s) and its r— derivative are 
bounded, while G k (s,r)F k (s)/F k (r) and its r-derivative are bounded by c^k 1 / 2 ^ 1 / 2 
and c*fc'/ 2+1 / 2 respectively for any r (cf. Lemma l43|). Further |s(s)/s(r)| ^ (1 — Si) 
and from (lllll). we have 



2fc-2+T 



fcl/2+5/2 ' 



and the lemma follows. 



I 



Remark 45. Since for r G (0, e], the bound in Lemma \J4\ on A\ is 0(k 
will see later that A k is dominated by A® ( defined in (llty) ) as k — > oo . 



we 



30 



O. COSTIN, J. L. LEBOWITZ AND S. TANVEER 



Lemma 46. Define GoiCiV) = fi m fc-s-oo G(C, v) an d Hq(C) = li m fc->oo H(C)> where 
C, <C A: 1 / 2 as /c — > oo. Then, 



(124) 



(125) 



"e-^fc(C,?)f@*7 
C ^o(C) 



c 



#o(C) 



#o(C) 



Proof. Using (|120p and (j 121 [) and the behavior of Bessel functions for small argu- 
ment, [J, it follows that for £, 77 <C A: 1 / 2 we have 

J+i/— i _ W+i„-i 
(126) g (C, fj) - lim g(C, ry) = 77 \ \ ? 



and H (Q = limfe-j.oo H(0 — J f C^e^z+iMO- Now, using the modified Bessel 
function equation, it is easily verified that /(C) — e~^H (() satisfies 

1(1 + 1) 



f" 



c 2 



/ = / 



with /(C) ~ e * as C — > oo. Using variation of parameters to invert the left hand 
side of the above equation, and using the boundary conditions at oo we obtain 



/(C) 



Go((,r))mdr) 



Dividing through by /(C) > the first identity in the Lemma follows. By differentiating 
the first identity with respect to Cj an d using the first identity in the resulting 
expression, we obtain the second identity. I 



Lemma 47. For any r £ (0, 1), 



(127) 



A k [l](r)-1 



1 Q( s ) rnk =j^G k (r,s)ds-l 
m k (r) 



For jr < r < | we get 



(128) ^MW, 
while for any r €E [0, 5], 
(129) 



1 

n(s) 



d 

dr 



d_ 
dr 



r -—G k (r,s)ds 

m k -i(s) 



k 3 r 



G k (r,s) 



ds ^ k 



m k (r) 

Proof. Recalling the definition (j9"3")l , it follows from and Lemma I4T71 that 



(130) 



C k m k - ilm k _ 1 = m k 



where j k (r) = 0(1) as k — > +00 for any r G [0, 1]. We can check from (f3"6")l that 
m k (l) — 0, m' k (l) — for k ^ 1. From ()130[) . inversion of yields 

r-l 



(131) 



m k (r) 



G k (r,s) < n(s)m k _ 1 (s) 



3k{s) 

5 (s) 



m k (s) } ds 
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Therefore, 



(132) 



1 Gk{rjS) mm p M ds = 1 _ r 1 Gk{rjS) j^M ds 



m k {r) 



s(s)m k (r) 



First, we choose Si so that 1 — Si = (5 + i)logfe/ (4fe + 2r). We then define 
5 so that (1 — Si)s(r) = s(r + S). It is clear that for large k we have S ~ 
(5/2 + l/2)s(r)logk/ ((2k + T)y/Q(r)). LemmaO and the fact that fc'/2-i/2(i_ 
<5i) 2fe+1 +7(2fc + 1 + r) < £ give 

(133) 



1 ^ / ~, jk(s)m k {s) 
G k (r, s)^^ — as 



r+<5 



s(s)m k (r) 
F k (r(st) 



1-5! 



2k+T 



Now, consider the contribution from J r+S ■ There are again two cases: (i) 1 r ^ 
fc- 2 / 3 and (ii) 0<r^ Ar 2 / 3 . 

In the first case, Taylor expanding G k (r,s) near s — r we get Gfc = (s — r) + 
0((s - r) 3 Q fc ) = - f) + 0(fc 4 / 3 (l - i) 3 , (1 - t) 2 ). Hence, 



(134) 



S(s)m k (r) 



< C*l|jfc|| 



t Ak+T -\l-t)dt < ~ 

K 2 



For the case (ii), we rewrite the integral in terms of £ = fear, to obtain 
(135) 



s(s)m k (r) 



^ L2 1 1 oo 



C+fcQ e-Q(,) + Q(C)g (Cr/) ||) d?? 



A- 2 



C+kaS 



by Lemma 051 Using (|132p and (|135l) . the first part follows. 

To prove (|128p , we note that if C3 is large and rk > C3, Taylor expansion gives 



(136) Ui{s,t) 



F k (r(st) 



y/Q(r(st))F k (r(s) 



G k (r(s),r(st)) 



U( S )(l-t) + 0((l-t) 2 ,k(l-tf, 



(1-tf (1-tf 



r 2 ' kr 2 



r l!s ./) = ,/>ki - n + o { (1 - *) 2 , k(i - tf, 



for / 4 = — s/fi(r(s)), while 

d_ 
ds 

From (I132P we note that 



(137, l M m . -^{(^ Ji|f.c 



'l-c5i 



2fe+T . 



3k{s)F k {s) 
drJr+g\t(r)J £{s)F k (r) 



G k (r,s)ds 
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We note further that 



(138) 



dr 



r+8 



2k+T 



jk{s)F k (s) 
mF k (r) 



G k (r, s)ds 



r+S \€(r) 



Or 

2k+T 



F k (s) 



£00 



Fk(r) 
jk{s)F k (s) 



G k (r,s) 



ds 



+ (2fc + r)^f/ ^ •>/ ; G k (r,s)ds 

/ A \ 2k+T 

, (ar + S)\ Mr + S)F k (r + 6l Gk ^ r + s 1 + Sl{r: 



S(r) 



Z(r + 5)F k (r) 



From the bounds in Lemmas l40l and l43l and the fact that £(s)/£(r) < (1 — 6%), we 
easily conclude that the contribution of f g in (|138l) to ^.A/.[l](r) is 0(l/k 2 ). 

Since Lemma l40l implies |jfe(r)| < c* and |ij.(r)| < + c*/(£x 2 ) for | > r > 
it follows from the local expansion of J7i(s, t) and its s-derivative in a neighborhood 
of t = 1 in the first integral in (|137[) that 



i r +5 Gfc(rS )MfW(£) ds 

*i r ' s(s)m k (r) 



fc 3 r 2 



jfe 5 



and (|T2"g)) follows. 

We now prove (|129|) . We first note that for r > fc~ 2 / 3 , s <E (r,r + 5), from 
pjf . d r G(r,s) = -1 at s = r and therefore, from (TL2T)]) , (fT2Tj) , it follows that for 
s - r = 0(/s _1 logfe) < /c~ 1/2 , <9 r G(r, s) — -1 < for s e (r, r + 5). The same is 
true for r € [0, fc~ 2 / 3 ] since in this regime, d r G k (r, s) ~ 9^Qq((, rj) (see (|126|l ). with 
C = r/(ak), rj = s/(ak). Therefore, from (|3Tj) and (f36|) . we get 



(139) 



d_ ( m k -i(r)G k (r, s) 
dr \ m k (r) 



(2k + r- 2) 



VW) F k (r) 



m fe _i(r)G fc (r, s) 



H(r) F k (r) 

- d r G k (r,s 



m k {r) 

m k -i(s) 
m k {r) 



Since the contributions to the integrals from J^ + g is 0{^), and the first term on 
the right on (|139|) is negative for large k, while the second is positive, it follows that 



(140) 



fi(s) 



d 



m k -i(r) 



dr \m k (r)G k (r,s) 



ds 



< \-rMm\ 

dr 



(2k + r- 2) 



Kir) 
Fk(r) 



\A k [l](r)\+0[~ ) _ < ■ k 
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for r G \Cik l] . For r G [0, C^k 1 ] , we note that since the contribution from 
S r +5 for i-Ak{l](r) is negligible, we have 



(141) ^[l](r)~to| 



C+fca5i 



-Q(»?)+Q(C) 



c 



oi \ ff(r/(l - l/k)) 



H(0 



ka- 



Ho(C) < 



#o(C)' 



it follows immediately from Lemma |4"6"1 that in this case , |^;^4fe[l](^)| < c*fc. Hence 
the inequality in (|140jl is valid for all r € [0, 1/2]. | 



Lemma 48. for any / e L°°[0, 1], 

(142) (a) for re [0,1], ||A/|| 



k 2 



11/11 



(143) 



(b) For r £ 



0, 



dr 



[A/](r) 



< c*fc||/|| 



Proof. Consider the expression for A k f from (|93j) . We break up the integral into 
J r r+A and jl +s , where 5 — C^fc -1 logfc, with C'2 large enough so that 



(1-S 1 ) 



From p6p and Lemma 1431 part (2), transforming the integration variable to t, it 
follows that 

(144) 



Q( S )^G,( M )/(# 

r+<5 Wfc(r) 



In f^ +S (we replace upper limit r + 5 by 1 if r + 5 > 1). Since 61 — O (k 1 log A;) 
and te(l- <5i, 1) then T fe (s,t) > and G fe (r, s) for r e [fc~ 2 / 3 , 1] Therefore, 

^• +a n(a)m fc _i(s) j 



(145) 



IIA/Hoo < 11/11 



m fe (r) 



-G fe (r,s) 



k 2 



From ()144|) we get we have 



nl ,m fe _i(s) 
S2(s) Gfe(r, s)as 



r+5 



rn fe (r) 



k 2 



Hence 
(146) 



n(«; 



rafc-i(s) 
m k (r) 



G k (r, s)ds 



Using Lemma H3 (|142[) (a) follows. For (b) we write 
d 



(147) 



dr 



m k (r) 



G k {r,s) 



mk-i(s) 
Wfc(r) 



f(s)ds 



By Lemma l47"l the quantity above is bounded by c*A;||/||oo. I 
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Lemma 49. For any f e [0, 1], 



||ftfc/||oo < pll/Hoo 



J^/KOIU ^||/||oo 

Proof. As before, we choose <$ = Cafc^logfc large C2 independent of fc. Using 
Lemma FEU it follows that 



(148) 



1 n(s)mk+l{s) G k (r, S )f(s)ds <C»(l-<Ji) 2k+a fc , / 2 - B / 2 ||/|| c<Sl|/|| 

r+5 ^fc(r) « 



(149) 



1 9 f f2(s)mfc+i(s) 



r+<5 



<9r [ m k (r 



G k (r, s)\f(s)d 



<cw(l-* 1 ) 2fc+a * , / 2 - 3 / 2 ||/|U<p||/|| 



Now, Lemma H21 implies 



(150) 



r+5 



n(s)m k+1 (s) 
m k ( r ) 



G k {r, s)f(s)ds 



c.J!4t [\ 2k+2+T dt < c * 11/11 



K JO 



A- 3 



(151) 



9 ffi(s)ra fe+ i(s) 



G k (r,s)U(s)ds 



k 2 



I 



Lemma 50. There exist ka and c*, independent of k, so that for k > ko, over the 
r -interval (0, 1), 



(152) 



\\hk\\oo < c* 



Proof. First we note that for k$ sufficiently large, ||/ifc ||oo exists since g k[) is contin- 
uous for r € [0, 1] and the expression for m k in (|36[) shows that l/m ko is bounded 
as well for sufficiently large ko since has no zeros in the region of interest. 

Define r k = H k h k+ i. Note that 

(153) h k = A k (Ak-ihk-2 + r k -i) +r k 

In k — fco inductive steps we get 

k— ko — 1 / m \ 

— m. h-k m + 1 

m=l y=i / 

We write this abstractly as 

(155) f) = fi°+91r, 
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where 

(156) \? k = A k A k -i..A ko+ ih ko ; 



k—kn — 1 / m 



— m^fe— m+1 

"i=l = 1 / 

and 9t is defined on the space S of sequences f) = {^a- }fcLfc +i m * ne norm 
(157) ||f)|| = sup \\hkWoo, 

fe^feo + l 

Lemmas 05] and |3U imply 
(158) 



k—ko — 1 I rn r 

»k Hoik ( p + m 

„ l= i I j—i 



1 



(fe - j + 1) 2 



(fc - m) 2 



<H|f)l 



where, if fco is large v < 1 is independent of fc. Thus, OT is contractive and there is 
a unique solution of ()155p in 5. | 

Lemma 51. _Fbr any r € [0, |] and for large enough k we have \\-£-,h k \\oo ^ c*fc. 
Proof. Since by Lemma [BUI fet. is bounded, Lemmas 1491 and 1481 imply 

\h' k (r)\ *S \^[A k h k ^](r)\ + \±[H k h k+1 ](r)\ < c*k 
ar ar 

I 

Lemma 52. For all k 1, /ifc(l) = 1. 

Proof. In case (i), a simple computation shows that 

92k g%r k \^o = i k h k (l); (g no - k t k m k h k ) 

(By the differential equation for h k , all derivatives exist.) Lemma 1371 with j = 2k 
gives 

Q s 2k I 

implying the result in case (i). In case (ii), using Lemma l38l a similar computation 
shows that 



ds 2k+ 

I 



Q2k+1 

i k = 2fc+1 U=o9n -k = i k h k (i) [g no -k ■= i k m k h k ) 



Definition 53. Let 

,Ac *\ = .«-2fc+l-r [%2k-2+T ( ^ 



(15?)) T k (s, S ) = S - 2k+1 - T I t 2k - 2+T s^-T k (s,t)dt, 



where T k (s,t) is defined in [95\ 
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Lemma 54. Let 5 = k UogA; and Sk{s) := ^ L t 2h 2 T k (s,t)dt. If C 2 is Zarge 
enough, s £ (0,5) and r(s) ^ fc _1 C*2 we have 



(160) T fe (s,s)=sS fc (s) 



O 



(1-s) 4 (1-s) 3 (1-s) 2 (1-s) (1-s) 3 (1-s) 5 



Proof. This simply follows by integrating (|103[) from t — 1 to s of T k and the fact 
thatf fe (s,l)=s5 fc (s). | 



4.13. Proof of Lemma 1311 First choose e! > 0. From Lemma |4"5I it follows that 
d 



ds 



[H k h 



fc+i. 



< p||/lfc+l||oo < p 



where we applied Lemma [5(11 Further, we note that 
1 d 



(161) 



(2/c + r)(2fc + T-l) 



Akhk-i(s) 



2fc+r-2 



We have 

(162) f t ^+r-2dT± 
Jo 



ds 

dt t 2k+ 



0e («,t)/i fc _i(Bt)dt + y < 2fe+T - 1 r fc (s,o^-i(st)^ 



(s,t)h k _ 1 (st)dt = h k -i(s)S k (s) 

S— (s,t) y h' k _ 1 (ss)ds = hk-i(s)Sk(s) 



T -2dT k 



h' k -i(ss) 



2k + r-2 s dTkr f)dt 
OS 



ds 



h k -i(s)S k (s) - h / k _ 1 (ss) S 2k - 1+T T k (s, s)ds = (2k + r- l)S k (s) 
Jo 

X / S^-^/lfe.iCss)^- / S 2fe - 1+T [f fc (B, S )-f fc (B,l)]/l' fc _ 1 (5s)ds 

Jo Jo 



Therefore, 



(163) 



(2fc + r)(2fc + T-l) 



[T fc (s,s)-T fe (s,s)+sS fc (s)]s ite 



2fc + T-l 



X h' k _ 1 (ss)ds+ (2k + r- l)S k (s) I * 

/o 



2fc+r-2 



h k ^ 1 (ss)ds 



We note that 



(2fc + r)(2fc + r - l)^(s) = | [^[1](b)] = O (~ ~ 
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and that (2k + r — 1) L s 2k+T 2 h k -i(ss)ds has a bound independent of fc. Com- 
bining (j!03p with Lemma [SH if k is large so that hex is large, then 

(164) T k ( S ,s)-[f k ( S ,s)-sS k (s)} = (l-s)+(- 1 ^ + ^ 



■{!-') 



-s - 



4 r 

Jj_ h 
12 3fcr 3 



(1- S ) 



O 



(1-s) 4 (1-s) 3 (1-s) 2 (1-s) (1-s) 3 (l-s) : 



kr 4 



k 2 r 4 



(1-s) 4 (1-s) 3 (1-s) 3 (1-s) 2 



kr 3 



kr 



From (HMD, h is clear that T fc (s, s) - T fe (s, s) + s,S fc (s) > if s € (1 - S, 1) and fcei 
is sufficiently large. Now, s/3/ (3fcr 3 ) (1 — s) 3 > exceeds any term following it in 
(|164p . except possibly when 1 — r, i.e. s is small. Thus, if we define 

(165) 



M k = sup \h' k {s)\ 

r(.e)e[ei,X] 



we get 
(166) 

\h' k (s)\ < {2k + r)(2fc + t- l)M fe _! 



„2fe+r-l 



l-<5i 



(1-s) 



4 



x [-I(l- s ) 2 + ?(l-s) 3 ]+ S g(l-s) 3 



rfs 



fc 2 fc 3 e 2 

When (1 — r) (and thus s) is small, we can replace the term sf[/ (12) (1 — s) 3 on the 
right side of the above equation simply by (1 — s) 3 , which is clearly bigger. From 
the fact that / 1 1 _ 4 s^-^-k' 1 (1 - s) 2 + (2/3)(l - s) 3 ]rfs = 0(fc~ 5 ), it follows that 



(167) 



2TTTT7 + p + ^ 2 J + p + ^ 2 



Let C3 be large enough and define fco(ei) = C%/e\, so that for fc ^ fco we have 



Then for fc ^ fco, 

(168) 
implying 

(169) M fe < 



2fc + r - 1 

2fc + r + 1 + ^fc 3 + fc 2 

1/2 



fc - 1 



1/2 



Af fc < 



1 



1/2 



fc 2 fc 3 e 2 

C* v -\ 1 ^ -s 1 

p72 ^3/2 + k l/ 2 2^ 7 -5/2 £ 2 



j=k 
I. 3 / 2 



: fc 1/2 fcV2^/ 2 fcl/2fc 3 / 2 e 2 

The result follows from the definition of M k and noting that last two terms in (|169p 
are CKcfc^V 1 / 2 ) 
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4.14. Proof of Lemma 1331 From Lemma l3i"1 and the definition of fco, it follows 
that 



-,3/2 

K(ei)\ < ■ " 



C4C3 




fc l/2 £ 3/2 

for fc ^ C^ei -1 = fcrj. Using /ife(l) = 1, it follows that for k > C^/r, 

1 r n 3/2 

Additionally, if HU afcr > — j = L c then \h k (r) - 1| < e. 

4.15. Proof of Lemma 1351 For £ 6 [0,L e ], using the a priori boundedness of hk 
in fc and Lemma EU we note that both hk(C) '■= ^fe( r (0) an( i (^fe)c are bounded 
independently of k. Hence the sequence {h k }k^2 is bounded and equicontinuous. 
By Ascoli-Arzela's theorem, there exists a subsequence hkj (C) converging to a con- 
tinuous function h. The first part of the result is proved. We first prove that 
IMC) - 1| < 4e. Now, from Lemma [551 

(170) |/ifc(C) — 1| < e for C € [L e , afc] for sufficiently large k 

Let /ifcj- be a subsequence that converges to h for £ g [0, L e ]. Let Cm, Cm be a 
minimum, and a maximum point of h on [0, L e ] and the corresponding minimum 
and maximum values are denoted by m and M respectively. Continuity at the 
endpoint ( — L e implies that M > 1 — e, rn < 1 + e. If both M — 1 — e < 
and m — 1 + e > 0, there is nothing to prove because in that case it is clear that 
|ft.(C) — 1| < 2e. Now, consider the possibility that (i): M > 1 + e. In a similar 
manner, we will also consider the possibility (ii): m < 1 — e. Consider (i) first. 
Since at the end point of the interval, h(L e ) < 1 + e, from continuity there exists 
an interval [a, b] C [Cm, L e ] of nonzero length for which 

(171) h(rj) < - (M + 1 + e) < M for 77 e [a, 6] 

For some L > L ei independent of k (to be determined shortly), we write 

(172) [Alf] (C)=(/ + ( J #(C,»7)/(»7(1 - fc" 1 ))^ 




[^/1(C) + K7KC) 



For fixed C and 77 we have 

(173) Urn K((, V )=K (t,r 1 ) = e-^^-g (<;,T 1 ) 



It is to be noted that for small enough e the inequality akr > L t always implies k > C3/V. 
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On our interval we have rj (. Thus Go ^ (see (|126[) ): Go can vanish only if 
i] — (. Furthermore, by (|171[) we have (m $ b]. We can then define 

3supr L i \h\ 

J = ~Fl Tj7 — > where K m = min K {CM,r)) > 

(b-a)K m ve[a,b] 

Note that Q(rj) ~ 77 for large fc and, aside from the exponential term, K is alge- 
braically bounded. We can thus choose L > L t large enough independently of k, 
so that 



(174) l^/lCOKeJ^H/IU,^,^] 

There is a subsequence of ^.fc j that converges uniformly on <E [0, L]; for simplicity, 
we will use the same notation h k j for the subsequence. It is clear that the limit 
is h(() if C € [0, L e ]. We keep the notation h for the limit on [0,L]. We note that 
(11701) implies 

(175) |MC) - IK e for ( G [L e ,L] 

Now choose a small 62 > 0. It is clear that in the interval [L e , L], h(() < 1 + e < M . 
For sufficiently large kj, using continuity of h(C), we have 

[A^HCm)} < f K(C M ,v)kv)dv + [ K((, t])h{rj)dr] + Me 2 

J'n£[(M,l]\la,b] Ja 

<M f K(( M , v)dv + \{M + 1 + e) / K(C M ,v)dr) + e 2 M 

Jv£[(M,L]\[a,b] Z J a 

= M [ K(Cm, V)dv ~ \{M - 1 - e) f K(( M , v)dv + Me 2 







< MAI°.[1](Cm) - ^^{M - 1 - e)K m + Me 2 

Since A kj [1] = A° k °.[l} + A% [1] + Al_ [1] (sec QUI and (JTT3J)) Lemmas S7J gU and 
(|174|) imply that for large fcj we have 



Hence, for large kj we have 



[A k j [l])(C M )<l + 1 +e 2 



e „ \ K 



a 



(176) [Al%](C M ) < M (l + - + 2e 2 ) - _J» (M - 1 - e)(6 

Now, there exists TV so that if j > iV, - Ml^^x] < e 2 and Aj = -4 fej+1 •••A J +i 
satisfies 

||Aj - /Hoc e 2 

while 

kj+i—kj — l m 

rj+x := B kj+1 + ^ X\_A kj+1 -i+iB kj+1 - m , 

m=l 1=1 

where Bi — %ihi + i, satisfies the estimate 
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Therefore, from 

hk ]+l = A,.4i /// +r J+ i 

it follows that 

hk ]+l (Cm) MCm) - £2 = M - e 2 
On the other hand, at £ = Cm we have 

(177) 

e _ , K, 



^jA kj h kj +r j+1 < (l+e 2 ) 
Thus, 

M - e 2 < (1 + e 2 ) 



M (1 + - + 2e 2 ) + e 2 ^ (M — 1 - e)(6 - a) 



+e 2 



M(l + - + 2e 2 ) + e 2 - ^(M- 1 - e)(6 - a) 



(2 



This is true for any e 2 , hence as e 2 4 0. Thus, 



M < 



M 



However, from the definition of J, this implies M — 1 — e ^ e. We note that 
for (ii), we repeat the above argument for —h, which has a maximum at £ m , to 
conclude that either (— m) — (—1 + e) < or (— m) — (—1 +e) = l — e — m^e. 
Therefore, 

l-2e<m^M<l + 2e 

implying that |ft — 1| < 4e. 

5. Appendix 

5.1. Short proof of the regularity of the unitary propagator. 

Theorem 5. Assume that Hi = H + V(x,t), where H is time independent and 
self-adjoint, and V(-,t) is in L 00 (W n ) for every t and is differentiate in time, with 
integrable derivative. Consider the Schrddinger problem 

(178) iip t = ip(x, 0) G D(H) 

Then there exists a strongly differ entiable unitary propagator on L 2 (K n ) U(t) so 
that ip(x,t) = U(t)ip € D(H) for all t and ip(x,t) solves [TTfy. 

Proof. We note that it enough to prove this property on a finite interval [0, e], since 
the problem can be restarted at t = e. Let y — ip — e - *?/V Then y satisfies the 
inhomogeneous Schrodinger equation 

(179) iyt = yae- 1 + Hy + Vy; y := iip + H^ + Vip 0) 2/(0) = 

We transform this equation into an integral equation, formally for now. Straight- 
forward calculations show that 

(180) ^ m y) t = e^e-tyo + e iHt Vy 
or (still formally) 

(181) ie lHt y = (J e (lR - 1)s ds^j y + j e lHs V{s)y(s)ds 

= (iH - l)" 1 ^*-* - 1)^ + /' e iHs V(s)y(s)ds 
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or, equivalently, 

(182) iy = (iH-l)- 1 {e- t -e- im )y + e~ tm I e ms V(s)y{s)ds 



o 

It is clear that (|182[) is contractive in the norm sup tg [ 0e ] || • || i 2( R 3) for small e, 
and has a unique solution. Clearly, the first term on the right side of (|182p is 
differentiable in time and the derivative is continuous since e~ lHt is; let uo denote 
this derivative. 

We now write a formal equation for u = y t . We have 

(183) iu = u +J^ e~ lHs V'(t-s) ^jf u(s')ds'\ ds+ e - lHs V{t-s)u(t-s)ds 

This equation is also contractive, and has a unique solution, in the same space. 
Thus both sides of (|183[) are integrable in time. By integration and appropriate 
changes of variables and order of integration, we see that f Q u(s)ds satisfies the 
same equation as y, which has a unique solution. Thus y — L u(s)ds is strongly 
differentiable. Since both y and e lHt y are strongly differentiable (the latter by 
inspection from (|18ip ). y £ D(H) for all t and is strongly differentiable. It is clear 
that ip £ D(H) and easy to check that it is differentiable and satisfies (|178p . | 

5.2. Laplace transform of the Schrodinger equation. We look more generally 
at equations of the form 

(184) iip t = Htp + V(t,x)ip 

where H is self-adjoint and time independent, and V(x,t) is bounded on R 3 and 
differentiable and bounded in t, and i/j(x, 0) £ D(H). The conditions on V can be 
relaxed. (For the purpose of this paper, H would be taken to be He-) 

Proposition 55. Under the assumptions above, the Laplace transform ip(p,-) of 
ip(t, •) exists for Rep > 0; it is in D(H) and satisfies 

(185) (p + iH)i> = V^o ~ iVip 

Proof. We take the unitary propagator of the time-independent problem, U = e~ lHt 
and apply U*(t) = U^it) to both sides of (fig!)) . Since (cf. flU]) C/" 1 is strongly 
differentiable, with derivative iU^ 1 !!, and ip is t— differentiable in L 2 , U" 1 ^ is 
differentiable and we get 

(186) (^"V)t = iU^HiP + U^tpt = -iU^Vip 

Since U~ 1 Vip is continuous in t, we can integrate both sides and get, after multi- 
plication by U and using the fact that t/ _1 (t) = U(—t), 

(187) i/j = Uip -iU [ U~ 1 Vi){s)ds = Uil)Q-i [ U(t - s)(Vip)(s)ds 

Jo Jo 

= Uipo ~iU* {Vtp) 

where * is the usual Laplace convolution. Taking the Laplace transform (which 
clearly exists) in (|187p and using standard functional calculus we get 

(188) ip= (p + iHy^jQ-iip + iH)- 1 ^ 
and thus ip is a D(H) solution of (|185p . | 
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Now, from eq. ([7]), it follows that y satisfies ©. Furthermore, using (|188[) and 
the fact that y° and flj are compactly supported, we see that y also satisfies 



(189) y(p,-) = < R X B y {p,-)-m X B 



where 9\q = (He — ip) 1 ■ 



^VLjyip-iju},-) 



5.3. Analyticity of (I — (£[ TO ) _1 in This is standard, and can be seen directly 
from analytic functional calculus. We provide a self-contained argument, for com- 
pleteness. We write £x to emphasize the X— dependence of £, and for simplicity 
of notation we drop the (I, m) subscript. We have 

(190) (I-£x 1 )- 1 -{I-£x>)~ 1 = (I-&x>)- 1 {£x 1 -Cx>){I-ex 1 )- 1 and 

{I - Cx')' 1 ^ + (€ Xl - € X ')(I - € Xl )~ 1 ] = (I-CxX 1 

We fix Xi and let X' -> X x . Since (I- £ Xl )~ 1 is bounded, then ||(C X i - £x>)(I ~ 
CxJ -1 !! as X' -> X 1 and 

(191) I + (£ Xl ~ £x-){I ~ CxJ -1 

is invertible when X\ and X' are close enough and [J+(£jti ~ £x')(I— Cxi) -1 ] -1 - ^ 
/ in operator norm as X' -> Xi. Thus 

(192) (7_ Cx ,)-i_^(/_ Cjfi )-i 

in operator norm, as Xi — > X'. Now diferentiability in X follows from (|190[) . 

5.4. Coulomb Green's function representation. The retarded Green's func- 
tions G = G+ is defined as the solution of the equation, 

(193) A G(x,x';k) = 5(x-x') 
in distributions, satisfying the radiation condition 

(194) G(x, x'; k) ~ F(0, <p)e lkl V" 1 ^; as r -> oo 
where 

(195) k=^/ip (Imfc > if Rep > 0), i/ = ^ 

Equivalently, G is the M 3 \{0} solution of (TTM|) with zero right hand side, satisfying 
(gSH and |x-x'|G(ar,a;';fc) -> (471-)- 1 as ar-a/ -> 0. 

Proposition 56. 

(196) ^oX B g= I G{x,x';k)g(x')dx' 

Jb 

Proof. The function 

(197) /:= / G{x,x'\k)g{x')dx' 



solves, as can be checked, the equation 

(198) A f = X B9 

with the radiation condition (|194[) . Such a solution is unique since the difference of 
two solutions satisfies the equation Aof — (with the radiation condition (11941) '). 
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Multiplying by G(x, x';k), integrating over a volume and passing to the limit where 
the volume approaches M 3 we see that / = 0. | 

Symmetries of the Coulomb potential — b/r allow for a closed form of G (cf. 26 - 
where the sign is chosen differently) in terms of Whittaker functions 2U and 9JT, 

(199) G(x; x'- k) = ^~ %V \, - W iv d-iH)m iv d-ikri) 



Anik\x~x'\ \d£ drj / 
where Im k > , 2kv = b and 

(200) f = M + W\ + \x - x'\ , J] = \x\ + \x'\ -\x- x'\ 

The Whittaker functions are defined in terms the Kummer functions M and U by 
the relations, see [1] Chapter 13, 

(201) m K>fl (z) = e~^z^ +fl M ^+ fj,- K,l + 2fi,zj , -7r<argz^7r 

%B K ^(z) = e~Sz5 +A1 [/ Q+/i-K,l + 2/i, zj , -7T < argz s$ it 

The following integral representation follows from [T] Chapter 13, for the values we 
are interested in, z% — ~ik^ 1 z-i = —ikrj 1 a = 1 — iu, b = 2 (a different "6" than the 
one in our Coulomb potential) 

™ / s e~i z zJ(z) , N e~i z zl(z) 

where / and J are as defined in (|5ip and the expression is valid in the regions where 
the integrals converge (in particular, |Im^| < 1). For other values of v of interest, 
the integrals can be replaced by appropriate contour integrals. For instance J would 
be replaced by 

(l - e - 2 ™) 1 I e zt r iv {l - t) iv dt 

where C is smooth simple curve encircling [0, 1], as it can be checked by calculating 
the jump across the cut of the integrand. It follows from these integral representa- 
tions that the Green's function is analytic at any (small) p, Rep ^ 0. Substituting 
QHD into COSO), we obtain 



5.5. Dependence of A in equation (1561) on Z, p. We now seek to determine 
the asymptotics of A in (1561) in the resolvent Xb^^Xb in terms of A = y/—ip and 

Z = cxp [wr&/(2A)] for X = (^/p, Z) E Df x D for sufficiently small e. 
Recall the expression A in (f55|). Note that since 

(203) a = ^X 2 - ic ~ e-"/ 4 c 1/2 (1 + 0(A 2 )) = a + A 2 ai + .. 

b 6e i7r / 4 

(204) Kl = 7T = ~ r" t 1 + °( A ')] = + + ' 
each of mi (a) and lOi(a) is analytic in A for small A, with the expansion 

(205) mi (a) = -Tt K l+ i/ 2 (2aa) ~ m 10 (a) +Xm 11 (a) + ... 

r 

(206) wi(a) = ~W Kl j +1/2 (2aa) ~ w lfi (a) + Xw ltl (a) + ... 
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The asymptotics in this case is also differentiable with respect to a and we get 
similar expressions as above for m'^a) and w[(a). It follows from the expression 
for /o in (|55p also possesses a regular series expansion in A: 

(207) /o(o) = /o,o(o) + A/ 0) i(o) + ... 

To simplify A as in (f56|) for small A, we now now consider the asymptotics of 102(a) 
and w 2 (a) for small A. 

5.6. Asymptotics of 102(a), w' 2 (a) for small A. Since 102(0) = ^2U K! ; + i/ 2 (2Aa), 
with k = b/(2X), it follows from formula (13.1.33) and analytic continuation to 
larger values of K of (13.2.5) of [T], p. 505 and the identity T(x)T(l — x) = n/ sin[7rx] 
that 

208) w 2 {a) = £ '- ^ '-H(2Xa; re, I) 

2ma 

where H(z; re, Z) = / e~ zt t l ~ K (l + t) l+K dt 
Jc 

where the contour C starts at ooe l °, circles around the origin once counter-clockwise 
to the right of t = —1 and goes to ooe l27r . In defining the integrand, we choose 
aigt £ [0, 2tt], arg(l + 1) € (— 7r, it] so that there is no branch cut on the real axis 
between —1 and 0. 

It follows from (|2"U5|) that 

209 w 2 (a = -A + - w 2 (a) + " -i?i(2Aa; re, Z) 

\ ay 27r« 

where H x (z;k,1) = [ e - zt t l - K+1 {l + t) l+K dt 
Jc 

We now seek to determine H(2Xa;b/(2X),l) and Hi (2Xa,b/(2X),l) asymptoti- 
cally for small A. For that purpose it is convenient to define 

(210) e 2 = 2A(^) 1/2 ,t = 6 2 ;, P( t; e 2 ) = -i log + + r, 

where we use the principal branch of log in defining P(r; e 2 ) above. Then, noting 
that in the definition of log r and log (r + e 2 ), argr G [0, 27r) and arg (r + a) S 
(— 7r, 7r], we have 

logT-log(r + e 2 ) = -log(l + — ) 

for r in the upper-half plane, while for r in the lower-half plane, we have 

£2' 



dr 



log r - log(r + e 2 ) = i2-n - log ( 1 
it is readily checked that 

b . 6' +1 / 2 



(211) tf(2Aa; A, |) = | T \ T + ^ e xp [-VabP(r; e 2 ) 

+ exp^-"^j J T l (T + e 2 ) l exp -VabP(T;e 2 ) 
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(212) g 1 (2Aa;A, = _^_ i ^ 



r i+1 (r + e 2 ) i exp \-Va~bP{ T - e 2 ) 



rir 



■ exp 



r i+1 (r + e 2 )' exp -Va6P(r; e 2 



c 2 



fir 



Here Ci is a contour in the upper-half complex r-plane from +00 to — e 2 along a 
steepest descent line, passing through the saddle point r Sl i = i(l + o(l)), where 
-P'(t S| i; 62) = 0. The contour C 2 is the steepest descent line in the lower-half r- 
plane from r = — e 2 to +00 through the saddle point r Si2 ,= — i(l + o(l)) where 
P'( T s.2: £2) — 0. We rewrite u> 2 and w 2 as 



(213) toa(o) = 

where 
(214) 

M 1 (C,e 2 ) = - 

7TZ 



(_ 1 )( + l e -Aa & /+l/2p (K _; ) r 



ZM^Vab, e 2 ) + Z _1 M 2 (2Va&, e 2 ) 



3 -CP(^ 2 



Ci 



V(r+e 2 ) z aY, M 2 (C,e 2 ) = - / e-^ T;e2 V(T+e 2 )W 



(215) 



w' 2 (a) = ( -A + z J w 2 (a)+ l 2 ;+i qA ; " ZM 3 (Vab) + Z^M^ab) 



a 



where 



(216) M 3 (C,e 2 ) = -[ e-' :P{T ^ ) T l+1 (T + e 2 ) l dT 
m Jc 1 

and M 4 (C,e 2 ) = -/ e- CP(T;£2 V +1 (r + e 2 )'rir 

It follows that, with e 2 = 2Xy/a/b, we have 

«4(fl) = _ A+ | (b 1 / 2 \ ( Z 2 A'h{Vab,e 2 )+Nh{Val>,e 2 ) 



W2{a) a ya 1 / 2 J yz 2 M 1 (Vab,e 2 ) + M 2 (Vab,e 2 ) J 

5.6.1. Analyticity in e 2 . 

Proposition 57. 77ie functions Mj(va6, •), i = 1, 4, are analytic near zero. 
Proof. We look at Mi , the others being similar. We can make a change of variable 
(218) q = P(r;e 2 ) -P(r M ;e 2 ), 

where the function g is real on the steepest descent contour and changes mono- 
tonically from 00 to 0, as we move from +00 to r = r S) i, and then increases 
monotonically again from to 00 as we move along the steepest descent path from 
t = t s 1 to r = — e 2 . We denote the two branches of the inverse function r(g) in 
dUBD by Ti(q) and r 2 (g). Noting that 
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we have 



(219) M l(C ,e 2 ) = e-^ (jf ^ ^ 



,z+i 



£2^2 



dq 



I r 1 2 + l + e 2 r 1 



It is easy to check that (t, — t Si ) 2 is analytic for small £2,, regular in g and nonzero 
at e 2 — for all g. 

Furthermore, the integrands in (I219[) are clearly bounded by an L 1 function 
uniformly in e 2 (see (|210l) and (|218|l ). ensuring e2-analyticity of the integrals. | 



(220) 
1 



Returning to the original variable r we get 



Mi(Vab,0) = — 

TTl TTl 



Cl.i 



exp 



i(2l + 1)0 - 2\/a6sin 
J21+1 feVab) - i \y 21+1 (2Vafo) + G21+1 (2Vab 



dd 



and 
(221) 
1 



M 2 (Vab,0) = — I e 

V 7TZ 7c 2 



-v^(-i+r) r 2i dr 



exp 



i(2Z + 1)6* - 2Va6sin 



(16 



— J21+1 (Vab) + 1 [Y 2l+1 (2^6) + G 2 i+1 (iVab 
where J2/+1 and Y2/+1 are the usual Bessel functions of order 21 4- 1 and 

/"OG 



(222) g 2 i+M = 



1 



(223) 



{exp[(2Z + l)t] + (-1) 2 ' +1 exp[-(2Z + } e 



sinh ((2Z + l)t) e 



J"... i = i\j 1 - -p - ~ = i + Series in e 2 
4 2 



„— v sinh £ 



-u sinh £ 



(224) 



Thus, asymptotically, to the leading order in A, we have with v = 2yab 1 
a w' 2 (a) 



b w 2 (a) 



(225) 



[Z 2 {J 2 i +2 (u) - iY 2i+2 (v) - iQ 2 i +2 {y)) + (J 2 i +2 (v) + iY 2 i +2 {v) + iQ 2 i +2 {v))] 
[Z 2 (J2i+i(v) - iY 2 i+i{v) - iG2i+i{v)) + {J 2 i+i(v) + iY 2l+1 (v) + iQ 2 i + i{v))} 

x (1 + 0(A)) 

The discussion on w' 2 {a)/w(a) shows that 

f (a)w' 2 {a) - fo(a)w 2 (a) 



A 



m' 1 (a)w 2 (a) — m\{a)w' 2 {a) 



is an analytic function of the extended parameter set X for X = (y/p, Z) G 
as long as the denominator for A is nonvanishing as A — > 0. 
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We can prove it is nonvanishing by simplifying the leading order expression in A 
for u> 2 (a)/u> 2 (a), defined as w' 2 (a)/w2 o( a ) under further assumption that a and c 
(as in the definition of /3) are sufficiently large. 

5.6.2. Further simplification for large a. For large a, there is additional simplifica- 
tion since 



(226) J a+ i foVab) ± iY 2i+1 (2Vab 



7r l/2 a l/4 5 l/4 



exp 



7r l/2 a l/4 6 l/4 



exp 



(227) J 2/+2 (2Vab) ± «r 2i+2 (2Va& 
and from Watson's Lemma, we get 

(228) G 2t+1 (2Va~b) = 0(l/o), fe+ 2 (2v^b) = 0(1 



-1/2- 



(Wa& + j 

(2Va7b-- 
V 4 



It follows that for large a, 

(229) ^oW.r^-Z 2 

where 

(230) n x 



w 2 ,o(a) a \Z 2 + rti 



= ie 4lVE " , Z = exp 



1 + 0(c 



iirb 
2A 



r 2 = i\ ab 



5.6.3. Nonvanishing of the denominator of A in i225\) . Now, defining 
(231) m = mi (a) , m! = mi (a) , / = f (a) , /' = f' a (a) 

We have to the leading order in A, for large a, 

" ni -Z 2 + 0(\) 
jh + Z 2 + 0{\a- 1 / 2 ) 



r-2 



(232) - A = 



-af 



r 2 m 



m - z 2 + o{\) _ 3 " 

th +Z 2 + 0(A,a-V2) ~ 

= / [4r 2 (m - Z 2 ) + 0(A)] - 4af [m + Z 2 + 0(A)] 
~ m [4r 2 (m - Z 2 ) + O(A)] - 4am' [m + Z 2 + 0(A)] 



The denominator of A is 



(233) 



D 



m\Z 2 Aa h 4r 2 + O(A) + m 4a 



4r 2 + 0(A) 



We note that 
(234) 



m = mi (a) = ;+1/2 (2aa) = e - QQ (2a)' + VM ( Z + 1 - 2Z + 2, 2aa ) 



(2A)' +1 a' 



r(i + i-^; 

and for large a in the fourth quadrant 



(235) 
(236) 



m' = m^ (a) ~ (2a) /+1 a' a 



r i + i 



[l + O (aa) 1 )] for a large 



[l + O^a)" 1 ] 



la, 



1 II 

ic^r c 1 exp 



as A -> 
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Therefore, D can be zero for large enough c (i.e. large /3) only if 

Z 2 (2V2ac 1 / 2 (l - i) [1 + O ((ca)" 1 )] + AiVb^j 

= -m (2\/2ac 1 / 2 (l - i) [1 + O ((ca)- 1 )] - Wfta) + 0(A) 
Taking the absolute square of both sides, we obtain, 

\Z\ 2 ([2aV2c] 2 (1 + Oic^a- 1 ) + [4Vab - 2aV2c (l + O^ 1 ^ 1 )] 2 ) 

= ([2aV2c] 2 (l + O^a" 1 ) + + 2a\f2~c(l + O^O] 2 ) + 0(A) 

This is impossible, since \Z\ < 1. This means that for large enough a and c (that 
is, j3 large), D cannot be zero. It means that the resolvent is well-defined as p = 
is approached from the closure of H. 

Note 58. Note that the denominator of A in i23ty) vanishes at points in the region 
\Z\ > 1, where, as a result, the resolvent 9\p has poles. From the relation between 
Z and p, it follows that p = is an accumulation point of a sequence of poles in 
the left half plane approaching zero tangentially to iR. 

5.7. Stationary phase analysis needed to calculate the ionization rate. 

We know that the solution y(is, x) is analytic in the extended parameter (yis, Z)~) , 
where 



(237) 
So, for X 

(238) 

Consider 
(239) 



is, z) e 



Z = exp [iirb/ (2y/s)] 



y(is,x) = J2s l / 2 F l (Z) 



1=0 



G(s) = Y,s l/2 F l (exp 



1=4 



iirb 



2^s 



It is clear that G(s) is a C 1 function of s in [—a, a]. Integration by parts gives 



(240) 

Now note that 
(241) 



G(s)e lst ds = 0{t- r ) 



Fi exp 



inb 
2Ts 



'2^1 



= E D i> 1 exp 

j>0 

with Djj decreasing exponentially with j, because of analyticity of Fi(Z) for \Z\ < 
1. For < I < 3, it follows there exists constants c and C independent of j so that 

3 

(242) ^2\D jtl \<Ce-^ 



1=0 



It follows that for large t, we have 



(243) 



IE E ^ 

i=i j=[Vt]+i 





' ibj ' 


/ exp 


.V*. 



e ist s l/2 ds\ < Cie- Cv/ * 
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e ist s l ' 2 ds 



< 



Co 



Further, for large t, 
(244) 
Therefore, 

(245) r (s' /2 ^(^)) e lst ds = ^ f" s l ' 2 F l {Z)e ist Us 



3 pa 

i=a J - a 



1=0 



o<;<3 



, , I s l/2 exp 



o<;<3 j=i 

We first evaluate the terms of the form 



i S St + J 



7T& 



ds + O 



(246) 

for large t, where 



J —a 

jnb 



dj = 
3 2 



The contribution from J° a is obviously small, at most 0(l/t), uniformly for all t, 
since the integrand vanishes exponentially as s — > 0~. So we only consider, for 

i<i< [Vt], 



(247) 



We have a point of stationary phase at s = sqj, where 



(248) 



so, 



2t 



2/3 



Note that soj <C 1 for t large since j is restricted to j < \fi. It is then convenient 
to rescale s — sojq, to obtain 



(249) s\y 2 exp [^(9 + 29- V2 



2-2/3 

g' /2 <ig, where 1/, = -^-i 1/3 



Using standard stationary phase arguments we obtain that, for large t, and hence 

large vj, 

(250) 



Jo 



l«n ,• / exp 
/o 



IV, 



q l ' 2 dq 



J 3 

For large i, the dominant contribution comes from the term with I = and so 



(251) 



[Vi] 



j y(is,x)e ist ds -^D jfi 



The sum over j is clearly convergent because of the exponential decay of -Dj.o! hence 
[Vt\ in the upper limit can be replaced by 00. From the definition of s ,j an d Vj, 
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it follows that 

(252) J y(is, x)e lst ds = O (V 5/6 ) 

At all other singular points, p = into, n G Z, the behavior is similar, and a 
similar calculation gives a e mwt t~ 5 / 6 contribution. Since Y g H, there is sufficient 
decay in n to ensure that the sum over all such contribution is convergent. 

5.8. Calculation of j^. Substituting the explicit expressions for mk(r) and m^-i) (r), 
it may be checked that in both cases, r = and r = 1, corresponding to (i) and 
(ii) respectively 

(253) jfe = fc 2 a 2 s(r)j| 2) + fcc^r)^ 1 ' +^ 0) , where 
(2) _ 417 ^ H(ak)H(C-C/k)\ , ff"(C) + 4v^)^'(0 



/ _ H(ak)H((-(/k) ^ H»(0 _ 1(1 + 1) 
a 2 s 2 \ H(a{k - 1))H(() J H (C) C 2 as(r) #«) 



(1) 2(l-2r)» / ff(afc)g(C-C/fc)) \ & 
' " " " * H(a(k-1))H(Q J a( 



2rVn us Q.' 



as 2v^ck 2an 



H(0 



(o) 5sO' 2 ujs 2 Q' sQ" (l + 2r) w 2 s 3 . 
J * = 4^-^^" S+ WQ- {UJn °- m)S 

where s(r) — y // Q(s)ds, £ = kar and jk := s[£femfe — Omfc-iJ/m/s. Recall that 
-ff(C) satisfies 

(2M) * SB % 1 "2fa? + 4fcan(0) H"? ^k) H ' 

where 

(255) a = 2-^ 

and that H(() has the following asymptotic behavior 

(256) ^(O-l + fcll + ^logC + O^ 1 ^,^ , (Cfc^oo, , C<M 

Now, we claim that for any r <G (0,1), \jk + I < Cfc~ 2 . In the regime 

r -C 1, we use Taylor expansion: 

(257) n = n(o) + ^'(o)^ + o (~\ , s - «(o) - VW)^ + o (~ 

and substitute r = (/(ka) in jf ) + fc" 1 ^^; we then use a = 2y/n(0)/s{0), (1254]) 
and the asymptotic behavior ()256|) to evaluate H(ak) and H(a(k — 1)) to find jj. + 
fc — Vfc 1 ^ ~ k~ 2 g(( t ) for some bounded differentiable function <?(£), with asymptotic 
behavior g(£) ~ const. /£ for large £. When r is not small, we use the asymptotic 
behavior (|256p to evaluate all terms involving the function H and to find the same 
inequality |j| 2) + k~ l j^ ] \ < Ck~ 2 . 

Therefore, jk(r) = 0(1) in all regimes. Further, it is easily checked that in the 
regime k >■ £ >■ 1, jfc(r) = 0(l,C _1 ) = 0(l/(fcr),l). Since the asymptotics 
is differentiable (since H satisfies a second order differential equation), it follows 
jl(r) = 0(k^ 1 r^ 2 , 1). When r is not small, using (|256l) . it is readily checked that 

f k = o(i). 
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5.9. Generalizations. In fact, the same asymptotic arguments hold more gener- 
ally if 

M 

V(t,x)= J2 e ij "%(r) 

j=-M 

with rij(r) satisfying the conditions we used for f2. We substitute for r = O(l) 



9n -k(r) 



T(2k/M + 1) 



exp 



M 



klogfo(r)+Y, kl ~ j/M M r ) 

3 = 1 



and calculate the error term R k as before. By requiring that the 0{k 2 ^ 2 ^ AI ) terms 
vanish for j = 0, .., M, we obtain (M + 1) first order differential equations for fj. 
To leading order 



Mr) 



\/fl- M (s)d£ 



The expressions for fj(r) for j ^ 1 are more complicated and involve arbitrary 
constants to be determined from the information for small k at r = 1. Again 
because of the presence of r~ 2 l(l + 1) in £k, the remainder is 0(r~ 2 ), which is 
0(k 2 ) when r = 0(k^ 1 ). We write 



(258) 



9n -k(r) ~ c* exp 



M 



fclog/oW + ^fc 1 -^ 1 ^-^ 



H(akr) 
T(2fc/M+ 1) 



Then, if C = 0(1), we find to leading order H(^) ~ -ffo(C) where 

#o _ 2 ^o — ^ " ' H ° = 

where now a = 2y/ fl^M (0) /s(0) and s(r) = J y / nZT/\7(s). As for Af = 1, we have 
to require H (() ~ 1 as C — > oo. This leads to 



; c C 1/2 ^ + i(C) 



#o(C) 



For nonzero g no -k, the constant multiple in (|258p is expected to be nonzero. On 
the other hand, the asymptotic behavior as ( J. 0, -ffo(C) ~ c*C~' implies that the 
behavior at r = of g nn -k/r is not acceptable unless every g„ vanishes identically. 

The analysis is likely to extend to systems with He replaced by 

Hw = -A - b/r + W{r) 

where b may be zero and W(r) = 0(r~ 1 ~ e ) for large r and is in L°°(R 3 ). Under 
these assumptions, W(r) does not participate in the asymptotics, to the orders 
relevant to the proofs. 

5.10. Further remarks on the asymptotics. 

Remark 59. A weaker statement than Theorem @] suffices to complete the proof of 
Theorem[TJ For instance, it suffices to show that for sufficiently large j, \Rk,j\ < 1, 
where 



fc .(r) =i^r l m kj (r)[l+R kj (r)] 



52 



O. COSTIN, J. L. LEBOWITZ AND S. TANVEER 



Remark 60. Stronger results than those in Proposition l36l hold. Noting that for 
any integer q ^ we have 



has the property limfc_>. 00 hk — 1. Indeed, this is in accordance with the heuristic 
arguments presented in <J5J3 While these results completely justify the formal 
asymptotics, they are not needed in the proofs and we omit the details. 
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